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Preface 


The  purpose  of  this  study  was  to  evaluate  one-dimensional  computer 
models  for  the  insulating  and  conducting  PRIZ.  The  PRIZ  optical  storage 
device  offers  great  promise  in  the  fields  of  optical  computing  and 
target  detection.  Both  the  Air  Force  Institute  of  Technology  (AFIT)  and 
the  Foreign  Technology  Division  (FTD)  have  sponsored  research  on  the 
PRIZ  since  the  early  1980's.  The  Soviet  Union,  however,  has  enjoyed,  or 
at  least  reported,  greater  success  in  experimentally  and  analytically 
studying  the  PRIZ  than  we.  This  thesis  is  one  of  several  produced  at 
AFIT  in  attempts  to  duplicate  reported  Soviet  results. 

Dr.  Theodore  Luke,  my  thesis  advisor,  has  coordinated  all  previous 
AFIT  investigations  into  the  PRIZ.  A  great  portion  of  my  time  was  spent 
simply  catching  up  with  the  wealth  of  information  he  has  accumulated 
over  the  years.  I  do  not  claim  to  understand  all  the  work  that  has  gone 
before  me,  but  I  credit  Dr.  Luke  with  helping  me  to  narrow  the  scope  of 
this  thesis  to  the  point  where  I  can  profess  some  expertise  in  my 
particular  subject  area.  I  especially  appreciate  his  patience  in 
guiding  a  computer-illiterate  helicopter  jockey  through  a  computational 
thesis.  I  also  appreciate  the  guidance  provided  by  Dr.  William  Bailey, 
one  of  my  thesis  committee  members,  who  often  helped  me  to  see  forest 
and  trees  at  the  same  time.  Finally,  I  thank  my  newl}wed  wife  Jennifer 
for  selectively  supporting  and  ignoring  me  at  the  appropriate  times. 

Joseph  W.  Cook,  III 


ii 


Table  of  Contents 


Page 

Preface .  ii 

List  of  Figures .  iv 

List  of  Tables .  viii 

List  of  Symbols .  ix 

Abstract .  xi 

I.  Introduction  .  1 

II.  Background .  9 

General  Description  .  9 

Theory  .  9 

Previous  Work .  17 

III.  Mathematical  and  Computer  Modelling  .  19 

Rate  Equations .  19 

Normalization  .  26 

Computer  Modelling  .  29 

IV.  Detailed  Analysis  of  Soviet  Journal  Articles  .  33 

Discussion  and  Interpretation  .  33 

Specific  Modelling  Goals  .  48 

V.  Results .  50 

Model  Validation .  50 

Reproduction  of  Low  Intensity  Exposure  Data  .  55 

Reproduction  of  High  Intensity  Exposure  Data  .  61 

Investigations  of  Reciprocity  .  70 

Investigations  of  Self  Erasure .  79 

VI.  Conclusions  and  Recommendations  .  88 

Appendix:  Landscape  Piets  .  90 

Bibliography  .  96 

Vita . 100 

iii 


List  of  Figures 


Figure  Page 

1.  PRIZ  Configuration  (Side  View)  .  10 

2 .  Energy  Band  Diagram  of  BSO  at  Room  Temperature  .  14 

3.  Typical  Dependencies  of  the  Intensity  R  of  the 
Readout  Image  on  the  Exposure  Energy  ^  for 
Various  Intensities  I'  of  the  Writing  Light: 

I'(l)  «  I'(2)  «  I'(3).  Taken  from  Figure  1, 

88  Paper  (6:374) .  37 

4.  Experimental  Dependencies  of  the  Maximum  Intensity 
^max  of  Readout  Image  and  the  Corresponding 
Exposure  Energy  ^^gx  on  the  Writing— Light 
Intensity  I'.  Taken  from  Figure  2,  88  Paper 

(6:375)  37 

5.  Net  Space  Charge  Density,  p,  versus  z  for  an  SI 
PRIZ  with  I'  =  50  pW/cm2,  iT  =  0.4  (1),  2  (2),  and 
4  pJ/cm^  (3).  Taken  from  Figure  3a,  88  Paper 

(6:375)  42 

6.  Net  Space  Charge  Density,  p,  versus  z  for  a 
conducting  PRIZ  with  n^  =  0.1,  1'  =50  pW/cm^, 

r  =  0.4  (1),  2  (2),  and  4  pJ/cm^  (3).  Taken  from 

Figure  3h ,  88  Paper  (6:375)  42 

7.  Net  Space  Charge  Density,  p,  versus  z  for  a 
Conducting  PRIZ  with  ng  =  0.1,  I'  =  500  mW/cm^, 

^  =4  (1),  8  (2),  12  (3),  and  24  pJ/cm^  (4). 

Taken  from  Figure  4,  88  Paper  (6:375) .  44 

8.  Comparison  of  Net  Space  Charge  Densities,  p,  for 

ng  =  0 . 1 ,  =4  pJ/cm^,  I'  =  50  pW/cm^  (1)  and 

I'  =  500  mW/cm^  (2).  Taken  from  Figures  3b  and  4, 

88  Paper  (6:375) .  44 

9.  Net  Charge  Accumulated  vs.  Time,  I'  =  50  pW/cm^, 

ng  =  0 ,  0.1,  Master  Model  One,  Parameter  Set  One  ...  54 

10.  Cathode  Electric  Field  vs.  Time,  I'  =  50  pW/cm^, 

ng  =  0,  0.1,  Master  Model  One,  Parameter  Set  One  ....  54 

11.  Photocurrent  vs.  Time,  I'  =  50  pW/cm^,  ng  =  0, 

0.1,  Master  Model  One,  Parameter  Set  One .  56 

iv 


12. 


Net  Charge  Density  vs.  Time,  Bin  1, 

I'  =  50  ^W/cm^,  ns  =  0 ,  0.1,  Master  Model  One, 

Parameter  Set  One .  56 

13.  Net  Charge  Density  vs.  Time,  Bin  2, 

I'  =  50  ns  -  0,  0.1,  Master  Model  One, 

Parameter  Set  One .  57 

14.  Net  Charge  Density  Distribution,  I'  -  50  pW/cm^, 

ns  =  0,  r'  =  0.4  (1),  2  (2),  and  4  (3)  .  58 

15.  Net  Charge  Density  Distribution,  1'  -  50  /iW/cm^, 

ns  =  0.1,  r*  =0.4  (1),  2  (2),  and  4  fxJ/cm^  (3)  ....  58 

16.  Effects  of  Absorption  Coefficient  and  Average 
Electron  Drift  Time  on  Net  Space  Charge  Density, 

I'  =50  /iW/cm^,  ns  =  0,  ^  -  4  /iJ/cm^  .  60 

17.  Comparison  of  Russian  Low  Intensity  Exposure  Data 
with  Our  Data  Scaled  up  by  a  Factor  of  26, 

I'  =  50  /iW/cm2,  ns  =  0,  ^  =  0.4  (1),  2  (2),  and 

4^J/cm2  (3)  (6:375) .  60 

18.  Comparison  of  Russian  Low  Intensity  Exposure  Data 
with  Our  Data  Scaled  up  by  a  Factor  of  26, 

1'  =  50  /iW/cm^,  ns  =  0.1,  r"  =  0.4  (1),  2  (2),  and 

4;iJ/cm2  (3)  (6:375) .  61 

19.  Net  Charge  Accumulated  vs.  Time,  I'  -  500  mW/cm^, 

Ug  =  0,  0.1,  Master  Model  One,  Parameter  Set  One  ....  63 

20.  Cathode  Electric  Field  vs.  Time,  I'  -  500  mW/cm^, 

ng  =  0,  0.1,  Master  Model  One,  Parameter  Set  One  ....  63 

21.  Photocurrent  vs.  Time,  I'  =  500  mW/cm^, 

ng  =  0 ,  0.1,  Master  Model  One,  Parameter  Set  One  ...  64 

22.  Net  Charge  Density  vs.  Time,  Bin  1, 

I'  =  500  mW/cm^,  ng  =  0,  0.1,  Master  Model  One, 

Parameter  Set  One .  64 

23.  Net  Charge  Density  vs.  Time,  Bin  2, 

I'  =  500  mW/cm^,  ng  =  0,  0.1,  Master  Model  One, 

Parameter  Set  One .  65 

24.  Net  Charge  Density  Distribution,  I'  -  500  mW/cm^, 
ng  =  0.1,  r  =  4  (1),  8  (2),  12  (3), 

and  24  (4)  nJ/cm^  .  .  .  ’. .  66 


V 


25. 


Comparison  of  Russian  High  Intensity  Exposure  Data 


with  Our  Data  Scaled  up  by  a  Factor  of  7, 

I'  =  500  mW/cm2,  n^  -  0.1,  r"  -  4  (1) .  8  (2)  , 

12  (3),  and  24  /iJ/cra^  (4)  (6:375) .  66 

26 .  Comparison  of  High  Intensity  Exposure  Data  for  an 
SI  and  a  Conducting  PRIZ,  I'  —  500  mW/cm^, 

ns  =  0,  0.1.  r  =  4  (1),  8  (2).  12  (3).  and 

24  ;iJ/cm^  (4)  .  68 

27 .  Comparison  of  Low  Intensity  Exposure  Data  for  an 
SI  and  a  Conducting  PRIZ,  I'  -  30  /iW/cm^, 

ns  =  0,  0.1,  ^  =  4  and  40  /iJ/cm^  .  68 

28.  Comparison  of  Medium  Intensity  Exposure  Data  for 
an  SI  and  a  Conducting  PRIZ,  I'  —  3  mW/cm^, 

ns  =  0,  0.1,  ^  ‘=4  and  40  /iJ/cm^  .  69 

29.  Comparison  of  High  Intensity  Exposure  Data  for  an 
SI  and  a  Conducting  PRIZ,  I'  =  300  mW/cm^, 

ns  =  0,  0.1,^  =  4  and  40  ^J/cm^  .  69 

30.  Relaxation  of  Net  Positive  Space  Charge  Following 
Turn-off  of  Write  Light,  I'  =  500  mW/cm^,  ng  =  0, 

r-  =  24  M-J/cm^ .  72 

31.  Relaxation  of  Net  Positive  Spate  Charge  Following 
Turn-off  of  Write  Light,  I'  =  500  mW/cm^, 

ns  =  0.1.  r  =24  pJ/cm^ .  72 

32.  Relaxation  Distribution  of  Net  Space  Charge  for 
High  (500  mW/cm^)  and  Low  (50  /iW/cm^)  Intensity 

Exposures,  ns  =  0 . 1 ,  =4  /iJ/cm^  .  73 

33.  Distribution  of  Net  Space  Charge  in  an  SI  PRIZ 
with  no  Delay  in  Recording  Response,  ^  =  4  (a), 

12  (b)  ,  20  (c)  ,  and  40  /iJ/cm^  (d)  .  75 

34.  Distribution  of  Net  Space  Charge  in  an  SI  PRIZ 
with  a  One  Time  Constant  Delay  in  Recording 
Response,  ^  =4  (a),  12  (b) ,  20  (c) ,  and 

40  /iJ/cm^  (d)  .  76 

35.  Distribution  of  Net  Space  Charge  in  a  Conducting 
PRIZ  with  no  Delay  in  Recording  Response, 

ns  =  0.1,  r  =4  (a),  12  (b) ,  20  (c) .  and 

40  ^J/cm^  (d)  .  77 


vi 


36.  Distribution  of  Net  Space  Charge  in  a  Conducting 
PRIZ  with  a  One  Time  Constant  Delay  in  Recording 
Response,  n^  =  0.1,  ^  =  k  (a),  12  (b) ,  20  (c) , 

and  40  /iJ/cm^  (d)  .  78 

37.  Net  Charge  Density  vs.  Time  in  Bins  1  (1),  2  (2), 

3,  (3),  and  4  (4),  I'  -  300  mW/cm^,  ~  0.1, 

^  =  0  to  4.9  mJ/cm^ .  82 

38.  Net  Charge  Density  vs.  Time  in  Bins  1  (1),  2  (2), 

3  (3).  and  4  (4),  I'  =  300  mW/cm^,  n^  -  10, 

^  =  0  to  4.9  mJ/cm^  .  82 

39.  Net  Charge  Density  vs.  Time  in  Bins  1  (1),  2  (2), 

3  (3),  and  4  (4),  I'  =  30  /iW/cm^,  nj  =  0.1, 

=  0  to  80  /iJ/cm^,  Master  Model  One  .  84 

40.  Net  Charge  Density  vs.  Time  in  Bins  1  (1),  2  (2), 

3  (3),  and  4  (4),  I'  =  30  /iW/cm^,  ng  =  0.1, 

^  =  0  to  80  ^J/cm^,  Master  Model  Two  .  84 

41.  Density  of  Ionized  Donors  vs.  Time  and  Distance 
for  SI,  Conducting,  and  DI  PRIZ,  I'  =  50  /xW/cm^, 

ns  =  0,  0.1,  ^  =  0to4  fiJ/cm^  .  91 

42.  Density  of  Ionized  Traps  vs.  Time  and  Distance,  SI 

PRIZ,  r  =  50  MW/cm2,  ns  =  0,  =  0  to  4  /iJ/cm^  ....  92 

43.  Density  of  Ionized  Traps  vs.  Time  and  Distance, 

Conducting  PRIZ,  I'  =  50  jjW/cm^,  =  0.1, 

^  =  0  to  4  ^J/cm^ .  93 

44.  Density  of  Free  Electrons  vs.  Time  and  Distance, 

SI  PRIZ,  I'  -  50  ng  =  0,  ^  =  0  to  4 

AiJ/cm^ .  94 

45.  Density  of  Free  Electrons  vs.  Time  and  Distance, 

Conducting  PRIZ,  I'  =  50  //W/cm^,  ng  =  0.1, 

=  0  to  4  ^J/cm^ .  95 


VI  1 


List  of  Tables 


Table  Page 

1.  Normalization  Factors  and  Normalized  Variables  .  27 

2.  Parameter  Sets . 51 


viii 


List  of  Symbols 


Roman 

Symbol  Definition 


D' 

do' 

r 

e 

E' 

Eo(t)' 

Ej(t)' 

g' 

G' 

H' 

h 

I' 

Jc' 

Jd' 
ddi  f  f 
Jt' 

n_' 

no' 

nbin 

Nd' 

He' 

neo(t) ' 

ns' 

Nt' 

t' 

to' 

texp 

H 

Vo' 


Electric  displacement  (C/cm^) 

Device  thickness  (cm) 

Total  exposure  energy  (J/cm^) 

Elementary  charge  (C) 

Electric  field  intensity  (V/cm) 

Electric  field  at  cathode  at  time  t  (V/cm) 
Electric  field  at  end  of  bin  j  at  time  t  (V/cm) 
Generation  coefficient  (cm'^*sec"^) 

Generation  rate  (sec'^) 

Magnetic  field  intensity  (Amp/cm) 

Planck's  constant  (J*sec) 

Write  light  intensity  (Watts/cm^) 

Conduction  current  density  (Amps/cm^) 
Displacement  current  density  (Amps/cm^) 
Diffusion  current  density  (Amps/cm^) 

Total  current  density  (Amps/cm^) 

Density  of  ionized  donors  (cm"^) 

Density  of  ionized,  or  filled,  traps  (cm"^) 
Concentration  normalization  factor  (cm”^) 

Number  of  discrete  and  equal  slices  of  dg' 
Density  of  donors  (cm“^) 

Density  of  free  electrons  (cm“^) 

Electron  density  at  the  cathode  at  time  t 
Injection  coefficient  of  cathode  (cm“^) 

Density  of  traps  (cm“^) 

Time  (sec) 

Initial  transient  time,  or  time  constant  (sec) 
Exposure  time  (sec) 

PRIZ  response 

Applied  external  voltage  (Volts) 

Longitudinal  distance  or  crystal  depth  (cm) 


Greek 

Symbol  Definition 


a' 


e 


Cd' 

Ct' 


Quantum  efficiency 
Absorption  coefficient  (cm~^) 

Permittivity  (F/m) 

Permittivity  of  free  space  (F/m) 

Relative  permittivity 

Donor  recombination  coefficient  (sec/cm^) 
Trap  recombination  coefficient  (sec/cm^) 
Mobility  of  electrons  (cm^/V»sec) 


IX 


p'  Charge  density  (C/cm^) 

Average  electron  drift  time  (sec) 

T’t'  Average  ionized  trap  lifetime  (sec) 

u’  Write  light  frequency  (Hz) 

The  prime  notation  (')  indicates  dimensional  variables  that  will  be 
normalized,  or  used  in  the  normalization  process.  The  normalization 
process  is  contained  in  Table  1,  page  27.  Dimensional  units  such  as  e, 
h,  e,  and  p  that  represent  well  known  physical  constants  or  specific 
material  parameters  do  not  carry  the  prime  notation.  Vector  quantities 
are  denoted  by  bars  as  in  E. 


X 


AFIT/GEP/ENG/91D-2 


Abstract 

This  study  attempted  to  analytically  reproduce  experimental  and 
analytical  data,  on  the  single-insulated  and  conducting  PRIZ,  presented 
in  1987  and  1988  by  the  Soviet  scientists  Bliznetsov  et  al .  Our  one¬ 
dimensional  analytical  models,  derived  from  earlier  Soviet  reports,  have 
consistently  reproduced  previously  reported  data  on  the  PRIZ.  In  this 
study,  the  same  codes  were  successfully  used  to  analyze  Che  effects  of 
injection  at  high  and  low  intensity  exposures  (30  /iW/cm^  to  500  mW/cra^)  . 
Although  the  codes  qualitatively  duplicated  the  nature  of  the  space 
charge  distribution  in  the  single-insulated  and  conducting  PRIZ,  our 
quantitative  results  were  significantly  different  from  the  Soviets' . 
Furthermore,  the  one-dimensional  codes  proved  insufficient  to 
effectively  duplicate  experimental  data  on  reciprocity,  self  erasure, 
and  regeneration.  Specifically,  we  were  unable  to  analytically 
reproduce  Soviet  experimental  data  establishing  a  reciprocity  region  for 
low  exposures  (<  20  /iJ/cm^)  in  the  conducting  PRIZ.  Since  the  basis  for 
PRIZ  operation,  the  linear  electrooptic  effect,  is  inherently  a  multi¬ 
dimensional  phenomenon,  we  were  not  surprised  to  find  some  limitations 
to  our  one-dimensional  models.  Our  results  should  serve  as  a  basis  for 
the  development  of  more  advanced  two-dimensional  analytical  models  of 
the  PRIZ. 
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ANALYSIS  OF  ONE-DIMENSIONAL  MODELS  OF  THE  INSULATING 
AND  CONDUCTING  PRIZ 

I .  Introduction 


The  PRIZ  is  an  optical  storage  device  that  offers  great  promise  in 
optical  computing  and  rmag3  processing.  The  term  PRIZ  is  a  Soviet 
acronym  for  Preobrasovatel  Izobrazheniya  which  translates  to  "linage 
transformer"  (32:343).  This  study  will  investigate  computer  models  of 
various  configurations  of  the  PRIZ.  The  output  from  these  models  will 
be  compared  to  previously  reported  experimental  and  analytical  work,  and 
the  results  of  these  comparisons  will  be  used  to  further  our 
understanding  of  the  basic  physical  processes  governing  the  operation  of 
t  . 

The  PRIZ  consists  of  a  single  cr/stal  of  photcrefractive  material 
•tounted  between  two  transparent  electrodes  which  are  biased  at  several 
kilovolts.  A  writ'' ng-tlight  beam  impinges  on  the  crystal  through  one  of 
the  transparent  electrodes.  In  general,  the  write  light  is 
perpendicular  to  the  crystal  face  and  impinges  from  the  direction  of  the 
negative  electrode,  although  other  configurations  have  been  used.^  The 
write  light  ioniz^'s  donor  sites  within  the  crystal  releasing  electrons 
that  drift  in  the  external  field,  producing  a  photocurrent.  Some  of  the 
conduction  electrons  are  captured  by  traps  within  the  crystal  prior  to 
reaching  the  anode.  The  ionization  of  fixed  donors,  drift  of 
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photoelectrons ,  and  capture  of  electrons  by  fixed  traps  produces  a  space 
charge  that  varies  with  time  and  distance  along  the  longitudinal  axis  of 
the  crystal.  The  s^ace  charge  will  also  vary  along  the  transverse 
dimensions  of  the  crystal,  depending  upon  the  characteristics  of  the 
writing-light  intensity  cross  section  and  length  of  exposure.  The 
transverse  electric  fields  that  are  produced  by  this  nonuniform  charge 
distribution  serve  to  modulate  a  reading-light  beam  that  follows  the 
writing-light  through  the  crystal.  The  read  light  is  of  much  lower 
photon  energy  than  the  write  light,  and  therefore,  is  assumed  to  have  an 
insignificant  effect  on  the  space  charge  produced  by  the  write  light. 

The  formation  of  a  space  charge,  which  produces  electrooptic  modulation 
of  the  read  light,  is  the  physical  basis  of  the  PRIZ. 

The  FRIZ  is  fabricated  in  three  configurations.  These  are  the 
double-insulating  (DI)  FRIZ,  the  single-insulating  (SI)  FRIZ,  and  the 
conducting  FRIZ.  In  the  DI  FRIZ,  both  electrodes  are  shielded  from  the 
crystal  by  transparent  insulators  that  prohibit  any  injection  or 
depletion  of  charge  within  the  crystal.  The  SI  FRIZ  has  an  insulator  at 
the  cathode  preventing  injection,  but  no  insulator  at  the  anode.  The 
conducting  FRIZ  has  no  insulators  and  therefore  allows  charge  injection 
and  extraction  into  and  from  the  crystal  through  the  external  circuit. 
These  differences  profoundly  influence  the  properties  of  the  FRIZ  in 
some  operating  modes  while  having  little  effect  in  other  modes.  To 
further  complicate  matters,  FRIZ  crystals  cut  from  different  bulk 
samples  can  exhibit  widely  varying  properties  (5:751). 
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As  mentioned  above,  the  PRIZ  has  several  operating  modes.  Though 
not  all-inclusive,  three  of  these  modes  are  image  framing,  dynamic  image 
selection  (DIS) ,  and  latent  imaging.  In  the  image  framing  mode,  which 
is  of  great  interest  in  optical  computing,  a  write  light  is  turned  on 
then  off.  After  a  suitable  time  for  image  development,  a  read  light  is 
turned  on,  and  the  resulting  modulated  image  is  observed.  The  read 
light  is  then  turned  off,  and  an  erase  beam  is  injected  into  the 
crystal,  while  the  electrodes  are  shorted,  to  redistribute  the  space 
charge  homogeneously  throughout  the  crystal.  The  complete  read-write- 
erase  cycle  is  called  a  frame.  As  reported  by  Casasent  et  al.  in 
"Applications  of  the  PRIZ  Light  Modulator",  write  times  can  be  as  short 
as  7  nsec  using  a  pulsed  laser  source.  Casasent  continues  with,  "In 
this  case,  photocarriers  were  generated  in  a  negligible  time,  but  the 
output  light  pattern  was  not  visible  until  1  ^sec  later,  and  it  peaked 
after  10  ^isec."  Casasent  also  reports  that  for  "more  conventional" 
exposure  times  of  1  msec,  the  output  image  is  "immediately  visible" 
(20:3847).  The  PRIZ  can  not  be  written  again  until  electrons  excited  by 
the  erase  beam  are  redistributed  and  relaxed.  Erasure  times  are  on  the 
order  of  1  msec,  and  one  millisecond  of  relaxation  time  is  generally 
adequate  before  commencing  a  new  frame  (20:3847). 

The  framing  mode  reveals  great  initial  insights  into  the  inner 
workings  of  the  PRIZ  and  the  time  intervals  involved.  In  experiments, 
Casasent  reports  repetition  rates  of  20  frames/sec  and  image  storage 
times  of  tens  of  seconds  to  one  minute  (20:3847).  In  a  survey  paper, 
the  Russian  scientist  Petrov  reports  repetition  rates  of  30  frames/sec 
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for  Bismuth  Silicate  (Bii2Si02o).  one  of  the  more  common  materials  for 
PRIZ  crystals  (32:345).  The  reader  who  is  interested  in  the  specific 
means  by  which  optical  computing  is  engineered  from  a  PRIZ  is  referred 
to  "The  Programming  System  PRIZ"  by  Mints  and  Tyugu  (27:359). 

A  second  mode  of  operation  is  dynamic  image  selection  (DIS) .  This 
mode  offers  great  potential  to  the  Air  Force  in  the  form  of  a  moving 
target  detector.  In  this  mode,  a  conducting  PRIZ  is  written  and  read 
simultaneously  and  continuously,  while  a  constant  voltage  bias  is 
applied  to  the  electrodes.  When  the  write  light  is  turned  on,  the 
modulation  of  the  read  light  peaks  and  then  decays  to  zero.  In  other 
words,  the  image  written  into  the  crystal  progressively  develops  and 
then  deteriorates  until  the  image  disappears .  Because  of  this  effect, 
images  of  stationary  objects,  which  present  a  constant  write  light 
intensity  to  the  crystal,  are  suppressed.  Images  that  introduce 
changing  intensities  to  the  crystal,  however,  are  enhanced.  This 
enhancement  can  increase  the  ratio  of  moving  image  intensity  to 
stationary  image  intensity  by  a  factor  of  20  (32:346).  The  PRIZ  in  DIS 
mode  will  also  regenerate  a  suppressed  stationary  image  when  the  write 
light  that  produced  the  image  is  turned  off.  The  underlying  physical 
basis  for  the  phenomena  observed  in  dynamic  image  selection  is  a  topic 
of  intense  investigation  and  will  be  studied  here  to  the  extent  one¬ 
dimensional  techniques  allow.  The  DIS  mode  of  operation  brought  the 
conducting  PRIZ  to  the  forefront  of  PRIZ  investigations. 

The  final  mode  of  operation  is  the  latent  imaging  mode.  In  this 
mode  a  latent  image  is  recorded  in  the  crystal  by  illuminating  the 
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crystal  while  the  voltage  bias  is  set  to  zero.  Although  the  resulting 
image  can  not  be  immediately  read,  subsequent  application  of  a  voltage 
bias  and  illumination  with  a  suitable  read  light  allows  the  latent  image 
to  emerge.  Since  this  mode  offers  little  additional  insights  to  the 
FRIZ  at  this  early  stage,  we  will  not  elaborate.  The  interested  reader 
is  referred  to  Petrov  (32:348). 

The  great  versatility  and  potential  of  the  FRIZ  should  be  obvious 
from  the  preceding  discussion.  Still  another  advantage  of  the  FRIZ  is 
its  inherent  low  noise  level  (32:343).  With  these  advantages  in  mind, 
however,  we  must  remember  that  the  FRIZ  is  a  very  complicated  device 
that  begs  intense  and  tedious  investigations.  One  very  significant 
problem  is  highlighted  well  by  Shlyagin,  et  al .  in  "Mechanisms  of 
Nonlinearity  in  a  FRIZ  Space-Time  Light  Modulator".  In  the  first 
paragraph  of  this  article,  he  states,  "Space-time  light  modulators  are 
used  in  optical  data  processing  systems  to  input  the  signals  to  be 
analyzed.  In  general,  it  is  desirable  for  the  control  and  output 
signals  of  the  modulator  to  be  linearly  related."  (37:68)  In  the  final 
paragraph,  he  confesses,  "We  have  thus  shown  that  in  general,  FRIZ 
space-time  light  modulators  transform  signals  nonlinearly . "  (37:72) 

Indeed,  studies  of  the  FRIZ  involve  nonlinearities  of  many  forms,  or 
more  specifically,  the  effect  of  nonlinearities  on  the  accurate  and 
consistent  production  of  optical  data.  For  example,  even  in  our 
relatively  simple  one-dimensional  models,  a  system  of  300  nonlinear 
differential  equations  is  used  to  predict  charge  distributions  within 
the  crystal.  Nonlinearities  are  therefore  at  the  forefront  of  FRIZ 
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investigation  and  development,  giving  the  device  its  great  promise,  but 
simultaneously  complicating  our  understanding  of  its  basic  operation. 

The  primary  goal  of  this  thesis  is  to  successfully  reproduce  and 
analyze  data  presented  by  the  Russian  scientists  A.  M.  Bliznetsov,  Yu. 

I.  Kuz'min,  and  A.  V.  Khomenko  in  their  paper,  "Investigation  of 
Reciprocity  in  the  FRIZ  Space-Time  Light  Modulator".  This  paper 
appeared  in  the  journal  Soviet  Physics  Technical  Physics  in  March  1988. 
Reference  will  also  be  made  to  self  erasure  results  reported  in  a 
previous  paper  published  in  the  same  journal.  This  earlier  paper  by  A. 
M.  Bliznetsov,  V.  V.  Bryksin,  L.  I.  Korovin,  S.  V.  Miridonov,  and  A.  V. 
Khomenko  entitled,  "Injection  Mechanism  of  Dynamic  Image  Selection  in 
PRIZ  Space-Time  Light  Modulators",  was  published  in  July  1987.  These 
papers  will  be  referred  to  as  the  1988  and  1987  papers,  respectively. 

Before  further  detailing  the  goals  of  this  thesis,  the  reader  is 
reminded  that,  although  the  PRIZ  derives  its  properties  from  the  multi¬ 
dimensional  aspects  of  the  electrooptic  effect,  one-dimensional 
modelling  is  the  first  step  in  achieving  a  basic  understanding  of  the 
PRIZ.  In  our  one-dimensional  model,  the  ultimate  effect  we  will  observe 
is  the  distribution  of  space  charge  versus  time  and  crystal  depth.  The 
space  charge  density  for  a  given  point  and  time  is  computed  by 
subtracting  the  density  of  conduction  electrons  (n^)  and  ionized,  or 
filled,  traps  (n_)  from  the  density  of  ionized  donors  (n+)  ,  and 
multiplying  the  result  by  the  basic  unit  of  charge,  e.  In  one 
dimension,  we  can  only  analyze  charge  distributions  produced  by  uniform 
illumination  of  the  PRIZ  face.  Such  illumination  produces  a  space 


6 


charge  that  can  only  vary  along  the  longitudinal  crystal  axis,  for  any 
given  time.  This  longitudinally  varying  space  charge  modifies  the 
longitudinal  electric  fields  within  the  crystal.  The  resulting 
transverse  fields,  however,  are  negligible  and  no  transverse 
electrooptic  effect  is  observed.  In  other  words,  the  read  light  will 
not  be  modulated  bv  a  PRIZ  under  uniform  illumination.  This  restriction 
severely  limits  our  ability  to  analyze  such  effects  as  reciprocity,  self 
erasure,  and  regeneration,  which  emerge  as  predominantly  multi¬ 
dimensional  effects  in  this  study.  Indeed,  some  well  documented 
experimental  effects  were  not  observed  at  all  in  our  one— dimensional 
mode  1 . 

Specifically,  this  study  will  use  well  known  electromagnetic 
formulas  and  a  differential  equation  solving  routine  to  predict  the 
charge  distributions  for  various  intensities  and  exposures  for  the  DI , 
SI,  and  conducting  PRIZ.  Using  this  scheme,  we  will  attempt  to 
duplicate  the  results  reported  in  Figures  Three  and  Four  of  the  1988 
paper,  Space  charge  data  will  be  manipulated  to  produce  information  on 
the  longitudinal  electric  field  and  current  densities  within  the 
crystal.  This  information  will  be  used  to  validate  the  authors' 
assertions  about  the  role  of  injection  currents  during  various 
illumination  intensities  and  exposures.  We  will  also  analyze  the 
authors'  assertions  about  reciprocity  in  the  PRIZ.  Finally,  we  will 
investigate  the  phenomenon  of  self  erasure  which  is  implied  in  Figure 
One  of  the  1988  paper  and  explicitly  studied  in  Figure  One  of  the  1987 
paper.  The  ultimate  goal,  of  course,  is  a  greater  understanding  of  the 
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basic  mechanisms  of  the  FRIZ.  This  paper  will  hopefully  exhaust  one¬ 
dimensional  efforts  to  unravel  the  FRIZ  and  provide  guidance  for  further 
multi-dimensional  models. 

The  following  sequence  will  be  followed  in  this  study.  First, 
background  information  will  be  presented  to  include  a  general 
description,  basic  theory,  and  previous  experimental  and  analytical 
work.  Next,  we  will  discuss  the  mathematical  and  computer  modelling 
process.  We  will  then  examine  the  1987  and  1988  papers  in  detail  and 
define  our  specific  modelling  goals  in  the  context  of  these  papers. 
Finally,  we  will  attempt  to  reproduce  the  experimental  and  analytical 
results  of  the  1988  paper.  Where  our  results  differ,  we  will  offer 
explanations,  and  where  our  one-dimensional  efforts  fall  short,  we  will 
highlight  multi-dimensional  studies  that  should  be  pursued. 
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II .  Background 


General  Description 

A  typical  PRIZ  consists  of  a  thin  slab  (approximately  one 
millimeter  or  less)  of  BSO  with  transverse  dimensions  of  10  mm  by  10  mm. 
The  two  large  flat  faces  of  the  crystal  are  coated  with  transparent 
insulators  in  the  case  of  the  DI  PRIZ,  while  only  the  flat  face  serving 
as  the  cathode  contact  is  insulated  for  an  SI  PRIZ.  The  conducting  PRIZ 
has  no  insulating  layers.  Transparent  electrodes  made  from  Platinum, 
In02,  or  another  conducting  material  are  deposited  over  the  insulating 
layers,  or  directly  onto  the  crystal.  These  electrodes  are  connected  to 
a  power  supply  that  can  provide  bias  voltages  of  several  kilovolts.  The 
write  light  can  be  laser  light  or  incoherent  light,  and  the  read  light 
is  generally  laser  light  of  a  much  smaller  photon  energy  than  the  write 
light.  The  write  light  normally  consists  of  an  intensity-modulated 
image  that  deposits  varying  amounts  of  energy  over  the  illuminated  face 
of  the  crystal.  In  our  one-dimensional  models,  however,  the  crystal  is 
uniformly  illuminated.  Appropriate  optical  devices  are  used  to  direct, 
focus,  and  analyze  the  write  and  read  beams.  Figure  1  shows  a 
simplified  schematic  of  a  PRIZ  device. 

Theory 

Basic  principles  from  the  fields  of  solid  state  physics, 
electromagnetism,  and  optics  are  used  to  understand  the  PRIZ.  When 
these  basics  are  combined  to  model  PRIZ  operation,  however,  the  combined 
effects  are  no  longer  basic,  and  the  macroscopic  results  are  often 
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Crystal 


Figure  1.  PRIZ  Configuration  (Side  View) 

difficult  to  predict.  This  section  will  provide  qualitative 
explanations  of  the  basic  physical  principles  governing  the  operation  of 
the  PRIZ.  The  mathematical  formulas  that  model  these  principles  will  be 
presented  in  section  III. 

The  first  effect  to  be  studied  is  photogeneration.  This  is  the 
process  by  which  donor  atoms  are  ionized  by  the  incident  write  light, 
releasing  electrons  into  the  conduction  band  of  the  PRIZ  crystal.  The 
primary  parameters  involved  are  donor  energy  level,  density  of  donors. 
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write  light  intensity  and  frequency  (or  energy) ,  quantum  efficiency,  and 
absorption  coefficient.  To  achieve  photogeneration,  the  energy  of  the 
write  light  must  be  at  least  as  great  as  the  donor  energy  level 
(expressed  in  units  below  the  conduction  band).  Given  this  condition, 
the  probability  that  a  single  donor  center  will  be  ionized,  if 
illuminated  by  a  photon,  is  expressed  by  the  quantum  efficiency.  The 
quantum  efficiency  used  throughout  this  report  is  one,  as  reported  by 
Sprague  for  BSO  illuminated  with  light  at  441  nm,  or  2.6  eV  photon 
energy  ( Sprague : 167 7 ) .  The  number  of  photoele^trons  generated,  by  light 
of  appropriate  energy,  will  generally  increase  with  increasing  values  of 
donor  density  and  write  beam  intensity.  However,  since  the  number  of 
donors  is  finite,  eventually  the  donors  will  be  depleted  and 
photogeneration  will  stop.  In  BSO  the  donors  are  Silicon  vacancy 
centers  located  2.6  eV  below  the  conduction  band.  The  absorption 
coefficient  characterizes  how  the  light  energy  is  absorbed  as  it  passes 
through  the  thickness  of  the  crystal.  Light  with  a  large  absorption 
coefficient  will  be  more  strongly  absorbed  in  the  initial  portions  of 
the  crystal . 

.More  complicated  additional  processes  exist  that  affect  the 
formation  of  a  space  charge  in  the  PRIZ  crystal.  For  example,  Peltier 
refers  to  luminescence  centers  at  1.3  eV  and  2.25  eV  in  BSO  that  can 
release  or  trap  conduction  electrons  (30:3684).  Furthermore,  we  know 
from  solid  state  physics  that  the  holes,  or  ionized  donor  centers,  can 
also  drift  under  the  influence  of  an  electric  field.  The  mobility  of 
holes  in  BSO  is  reported  to  be  so  low,  in  comparison  to  the  mobility  of 
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conduction  electrons,  as  to  be  inconsequential  in  our  calculations 
(14:1488).  Most  of  the  additional  processes  alluded  to  above  are  also 
relatively  minor  effects,  and  therefore,  are  not  specifically  addressed 
in  our  models  or  the  Russian  models. 

Once  the  photoelectrons  are  released  into  the  conduction  band, 
they  will  drift  in  the  external  field  produced  by  the  transparent 
electrodes.  The  rate  of  this  drift  will  depend  on  the  strength  of  the 
electric  field  and  the  electron  mobility,  a  material  property.  The 
drift  will  stop  when  the  electron  reaches  the  anode  or  when  it  is 
captured  by  a  trap  site  within  the  crystal.  In  BSO,  these  trap  sites 
can  be  shallow,  0.3  eV  to  0.6  eV,  or  deep,  1.3  eV  and  2.25  eV  (12:878, 
21:9).  The  deep  traps  correspond  to  the  luminescence  centers  mentioned 
earlier.  The  rate  of  trapping  will  depend  on  the  number  of  electrons  in 
the  conduction  band,  the  density  of  trap  sites,  and  the  trap 
recombination  coefficient.  The  trap  recombination  coefficient  is  a 
material  parameter  that  reflects  specific  quantum  characteristics  of  the 
traps.  Like  donors,  traps  are  finite  in  number  and  can  be  completely 
filled.  In  BSO,  the  average  time  between  photogeneration  and  trapping 
is  on  the  order  of  0.1  msec  (6:376). 

Since  the  donor  and  trap  sites  are  considered  fixed  within  the 
crystal,  when  donors  are  ionized  and  release  electrons  that  drift  to 
trap  sites,  the  result  is  a  net  space  charge  within  the  volume  of  the 
crystal.  At  any  given  point  within  the  crystal,  the  net  space  charge 
density  can  be  calculated  by  subtracting  the  charge  density  of  free 
electrons  and  negatively  ionized  traps  from  the  charge  density  of 
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positively  ionized  donors.  Neglecting  further  photogeneration  (the 
write  light  is  turned  off)  and  electron  injection  through  the  cathode; 
the  only  way  for  an  existing  space  charge  to  change  is  through 
deionization.  This  can  occur  through  recombination  of  ionized  donors 
with  conduction  electrons  or  through  detrapping,  in  which  a  trapped 
electron  is  thermally  or  optically  excited  back  into  the  conduction 
band.  These  two  processes  are  dependent  upon  the  donor  recombination 
coefficient  and  the  average  ionized  trap  lifetime,  respectively.  The 
donor  recombination  coefficient  is  a  quantum  mechanical  parameter  that, 
together  with  the  number  of  free  electrons  and  number  of  ionized  donors, 
determines  the  rate  at  which  the  ionized  donors  will  deionize  or 
recombine.  The  average  trap  lifetime  due  to  thermal  excitation  for  BSO 
is  very  large  at  standard  temperatures.  This  means  that  de trapping,  or 
deionization  of  traps,  is  a  minor  effect.  Detrapping  due  to  optical 
excitation  by  the  read  light  is  also  a  minor  effect  since  the  quantum 
efficiencies  involved  are  very  low.  This  is  in  spite  of  the  fact  that 
read  beam  energies  are  on  the  order  of  2  eV.  Figure  2  is  a  diagram  of 
the  energy  levels  in  BSO  at  room  temperature. 

The  previous  paragraph  detailed  the  individual  processes  by  which  a 
space  charge  is  formed  in  the  FRIZ  crystal.  The  formation  of  this  space 
charge,  which  varies  with  position  and  time,  is  the  basis  of  FRIZ 
operation.  It  is  this  space  charge  distribution  which  modifies  the 
longitudinal  electric  field  within  the  crystal.  In  the  more  general 
case  of  nonuniform  illumination,  the  space  charge  will  also  produce 
transverse  electric  fields  within  the  crystal  that  modulate  the  read 
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Conduction  Band 


Figure  2.  Energy  Band  Diagram  of  BSO  at  Room  Temperature 

beam  through  the  transverse  linear  electrooptic,  or  Pockels,  effect, 
which  will  be  discussed  later. 

The  macroscopic  result  of  space  charge  formation  in  time  is  a  very 
complicated  process,  however,  some  general  trends  can  be  presented.  The 
most  obvious  effect  is  the  formation  of  a  net  positive  space  charge 
region  in  the  body  of  the  crystal  next  to  the  negative,  illuminated 
electrode.  A  negative  space  charge  region  is  formed  near  the  anode  as 
the  photoelectrons  drift  under  the  influence  of  the  external  electric 
field  and  either  become  trapped  or  accumulate  at  the  anode.  Since  no 
insulator  is  present  on  the  anode  in  an  SI  or  conducting  FRIZ,  the 
electrons  that  drift  there  will  be  extracted.  The  resulting 
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longitudinal  charge  separation  modifies  the  electric  field  within  the 
crystal.  In  the  region  of  the  crystal  where  the  net  space  charge 
changes  from  positive  to  negative,  the  net  electric  field  will  be 
reduced.  As  a  result,  the  conduction  electrons  in  this  region  are 
slowed  and  experience  enhanced  trapping.  The  enhanced  trapping  further 
modifies  the  net  space  charge  and  the  net  electric  field  within  the 
crystal.  As  a  result,  the  region  of  enhanced  crapping  moves  toward  the 
cathode  as  photogeneration  continues.  As  the  net  positive  space  charge 
region  retreats  behind  this  region  of  enhanced  trapping,  the  slope  of 
the  net  space  charge  distribution  curve  at  the  ''rossover  point,  or 
bottleneck,  becomes  more  steep  and  r.aunced  (13:33,  32:336). 
Eventually,  the  net  positive  space  charge  is  confined  to  a  very  narrow 
region  adjacent  to  the  cathode.  The  process  is  modified  if  injection  of 
electrons  is  allowed  through  the  negative  cathode.  This  can  only  occur 
in  the  conducting  PRIZ,  and  the  specific  effects  on  the  space  charge  are 
a  primary  focus  of  this  report.  At  this  point,  we  can  assert,  however, 
that  the  injection  will  depend  on  the  material  properties  of  the 
electrode  and  the  net  electric  field  at  the  cathode.  Also,  in  the  case 
of  uniform  illumination,  we  can  expect  uniform  injection,  although  the 
injection  will  most  certainly  depend  on  time. 

The  PRIZ  is  actually  a  refinement  on  a  previous  optical  storage 
device  called  the  PROM.  PROM  is  an  acronym  for  Pockels  Readout  Optical 
Modulator  and  refers  to  an  optical  device  that  modulates  light  according 
to  the  longitudinal  linear  electrooptic  (Pockels)  effect.  This  effect 
occurs  when  a  longitudinal  electric  field  is  applied  to  a  suitable 
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crystal,  such  as  BSO,  in  an  appropriate  orientation.  The  electric  field 
modifies  the  crystal's  index  of  refraction,  producing  a  fast  and  a  slow 
axis  within  the  crystal  making  it  birefringent  (25:354,  33:816).  When 
light  traverses  a  birefringent  crystal,  phase  changes  are  induced. 

Since  the  originally  applied  electric  field  is  modified  within  the 
crystal  when  the  write  light  induces  a  space  charge,  the  birefringence, 
and  subsequent  net  phase  change,  is  very  complex. 

The  iJ  .3ar  electrooptic  (Pockels)  effect  described  above  is  a 
result  of  the  longitudinal  electric  field  within  the  PROM  crystal.  The 
PRIZ,  however,  modulates  light  according  to  the  transverse  electrooptic 
effect.  This  effect  is  due  to  the  transverse  electric  fields  produced 
by  the  space  charge  distribution.  These  transverse  fields  modulate  the 
read  light  as  it  passes  through  the  crystal.  The  longitudinal  fields  do 
not  contribute  a  net  modulation  to  the  light  (31:823,  32:330,  37:68). 

The  PRIZ  and  PROM  devices  are  similar  in  geometry  and  method  of 
operation,  however,  the  differences  in  modulation  mentioned  above  are 
produced  by  different  crystal  orientations.  As  stated  by  Gardner  and 
Luke,  "PRO.M  devices  use  a  (100)  crystal  orientation  in  which  the  [100] 
direction  is  normal  to  the  sandwich  structure.  PRIZ  devices  use  either 
Qll)  or  (110)  orientations  with  the  corresponding  direction  normal  to 
the  sandwich  structure."  (23:19)  The  sandwich  structure  refers  to  the 
crystal  at  the  center,  with  dielectric  insulators  on  either  side  (for  a 
DI  PRIZ),  and  transparent  electrodes  over  both  insulators.  Therefore, 
through  selective  orientation,  either  the  longitudinal  or  transverse 
electrooptic  effect  can  be  highlighted.  In  addition  to  BSO,  other 
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suitable  crystals  are  produced  from  Bismuth  Titanate  (Bii2Ti02o) ,  Bismuth 
Germanate  (Bii2Ge02o) ,  Lithium  Niobate  (LiNb03) ,  and  numerous  other 
compounds  (17:757).  Specific  modulation  of  a  suitable  crystal  will 
depend  on  individual  material  properties  such  as  electron  mobility, 
donor  energy  level  and  density,  density  of  traps,  crystal  symmetry,  and 
electrooptic  coefficients.  In  this  report,  we  will  deal  exclusively 
with  BSO  even  though  more  efficient  crystals  have  been  developed  and 
studied  (34:107).  It  is  interesting  to  note,  once  again,  that  uniform 
illumination  will  not  produce  significant  transverse  electric  fields 
within  the  PRIZ  crystal,  and  therefore,  the  transverse  electrooptic 
effect  will  not  be  observed  (20:3849).  This  highlights  a  limitation  of 
the  one-dimensional  study  offered  here. 

Previous  Work 

Much  of  the  analytical  and  experimental  work  on  the  PRIZ  has  been 
conducted  in  the  Soviet  Union.  Bryksin,  Petrov,  Bliznetsov,  Khomenko , 
Shlyagin,  and  Astratov  have  published  numerous  reports  on  the  PRIZ. 

-Much  of  AFIT's  work,  to  include  this  thesis,  have  concentrated  on 
experimentally  or  analytically  reproducing  Soviet  results.  These 
efforts  have  met  with  overall  success  as  reported  by  Shields,  Nilius, 
Cushing,  Gardner,  and  Nelson  (36:1,  29:1,  21:1,  22:1,  28:1).  Anderson 
provided  pioneering  work  on  specific  aspects  of  PRIZ  crystal  damage,  due 
to  high  internal  fields,  in  his  1986  AFIT  thesis  (1:1).  Luke  and 
Gardner  quantify  the  threshold  energy  density  for  damage  at  1.0  mJ/cm^ 
(23:21).  This  experimentally  derived  damage  threshold  should  be  kept  in 
mind  throughout  this  analytical  study,  as  we,  and  the  Russians,  are 
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often  dealing  with  exposure  energy  densities  well  above  these  values. 

Our  analytical  solutions  do  not  model  the  catastrophic  failures  that 
experimentally  occur. 

The  Russian  articles  leading  up  to  the  1987  and  1988  papers  studied 
in  this  report  show  varying  degrees  of  sophistication  and  accuracy. 

Some  papers  offer  gross  approximations  that  are  surprisingly  accurate  in 
some  operating  regimes.  For  example,  Bryksin  neglects  trapping 
altogether  in  a  1984  paper,  and  again  in  a  1986  paper  offering  an  "exact 
solution"  for  the  photoinduced  charge  in  photorefractive  crystals 
(10:2195,  9:80).  Similarly,  Bryksin  and  others  frequently  ignore  the 
absorption  coefficient,  and  assume  that  the  write  beam  is  absorbed 
totally  in  a  certain  crystal  thickness  near  the  cathode  (14:1488). 

Other  parameters  that  may  or  may  not  be  considered  are  depletion  of 
donors,  filling  of  traps,  thermal  ionization  of  traps,  recombination  of 
donors,  diffusion  current,  and  injection.  As  should  be  evident,  one  of 
our  major  problems  was  the  uncertainty  as  to  what  approximations  the 
Russians  used  in  the  particular  papers  we  studied.  This,  coupled  with 
imperfect  translation,  made  thorough  detective  work  as  invaluable  as 
mathematical  and  physical  rigor. 
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Ill .  Mathematical  and  Computer  Modelling 


As  previously  discussed,  the  mathematical  models,  and  subsequent 
computer  models,  used  depend  upon  the  approximations  made.  With  this  in 
mind,  we  have  developed  two  distinct  master  models  from  which  all  minor 
variations  are  applied.  The  first  master  model  is  comprehensive  and 
gleans  the  best  modelling  information  from  the  entire  body  of  previous 
Russian  and  AFIT  work  on  the  FRIZ.  Few  effects  are  neglected  in  this 
model  although  some  effects  are  numerically  insignificant  in  most 
operating  regimes.  For  example,  depletion  of  donors  is  modelled  even 
though  Bryksin  reports  that  depletion  of  donors  has  never  been  observed 
(18:2052).  Likewise,  thermal  ionization  of  traps  is  modelled  even 
though  the  characteristic  time  for  this  event  is  several  orders  of 
magnitude  greater  than  most  FRIZ  exposure  times. 

The  second  master  model,  or  clipped  model,  only  models  those 
effects  that  are  explicitly  mentioned  in  the  1988  paper  by  Bliznetsov, 
et  al .  These  are  trapping,  depletion  of  donors,  absorption  length 
modelling  for  the  write  beam,  and  cathode  field-dependent  injection 
(6:374).  All  models  used  in  this  study  will  be  based  on  one  of  the  two 
master  models  presented  above.  Minor  deviations  from  the  standard 
models  will  be  highlighted. 

Rate  Equations 

Each  of  the  m.aster  models  consists  of  three  basic  rate  equations. 
These  are  the  net  rate  of  production  of  ionized  donors,  the  net  rate  of 
production  of  ionized  traps,  and  the  time  rate  of  change  of  free 
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electrons  in  the  conduction  band  of  the  crystal.  When  these  equations, 
with  appropriate  initial  and  boundary  conditions,  are  applied  to  a 
specific  location  in  the  crystal  at  a  specific  time,  the  number  of 
ionized  donors,  ionized  traps,  and  free  electrons  can  be  calculated. 

The  determination  of  appropriate  boundary  conditions  is  nontrivial  since 
each  region's  boundary  conditions  are  dependent  on  the  other  regions' 
boundary  conditions.  With  this  information,  however,  the  net  space 
charge  density  can  be  calculated  along  with  the  net  electric  field  and 
current  density  for  that  particular  region  of  the  crystal.  The  same 
three  equations  with  different  initial  and  boundary  conditions  can  also 
be  used  to  compute  the  net  space  charge  density  in  the  other  regions  of 
the  crystal.  The  macroscopic  result  is  an  overall  view  of  the  net  space 
charge  density,  net  electric  field,  and  current  density  for  the  entire 
crystal.  The  space  charge  distribution  will  evolve  in  time  according  to 
the  rate  equations,  and  the  overall  result  can  be  used  to  predict  the 
longitudinal  and  transverse  fields  that  lead  to  the  longitudinal 
electrooptic  effect  (PROM)  and  the  transverse  electrooptic  effect 
(PRIZ) .  Accepting  the  rate  equations  as  valid  mathematical 
representations  of  the  physical  processes  within  the  PRIZ,  the  accuracy 
of  the  results  will  depend  on  the  accuracy  of  the  material  parameters 
entered  in  the  equations,  the  bin  size,  and  the  tolerance  of  the 
differential  equation  solving  technique  used.  The  bin  size  is  the 
fractional  slice  of  the  crystal  over  which  the  rate  equations  are 
applied.  Bin  size  will  be  discussed  in  detail  in  a  later  subsection. 
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The  first  rate  equation  is  the  net  rate  of  donor  ionization.  For 
master  model  one,  this  equation  models  exponential  absorption  of  light 
using  the  absorption  coefficient,  depletion  of  donors,  and  recombination 
of  donors.  The  corresponding  rate  equation  for  the  clipped  master  model 
does  not  include  the  recombination  of  donors  term.  Eqs  (1)  and  (2) 
below  give  the  net  donor  ionization  rate  equations  for  master  models  one 
and  two,  respectively: 


dn,  ' 
dt' 


n,')- 


(1) 


dn^' 

dt' 


n,') 


(2) 


where 


n+'  =  density  of  ionized  donors  (cm"^) 

t'  =  time  (sec) 

G'  =  generation  rate  (sec~^) 

N(j'  =  density  of  donors  (cm~^) 

ng'  =  density  of  free  electrons  (cm"^) 

Cjj'  =  donor  recombination  coefficient  (sec/cm^) 


The  prime  notation 


(')  for  the  variables  given  above  denotes  dimensional 


variables  that  will  undergo  normalization.  The  generation  rate,  G' , 


given  above  is  dependent  on  numerous  variables  as  described  in  section 


II.  Equation  (3)  on  the  following  page  gives  the  mathematical 
formulation  for  G' ,  where  we  have  introduced  a  new  parameter  called  the 


generation  coefficient,  g' : 
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a 


(3) 


G'  -  -  g'-^exp{-a 'z') 


hv  '  N. 


where 

I'  =  write  light  intensity  (Watts/cra^) 
r)  =  quantum  efficiency 
h  =  Planck's  constant  (J*sec) 
i/'  =  write  light  frequency  (Hz) 

a'  =  absorption  coefficient  (cm”^) 
z'  =  longitudinal  crystal  depth  (cm) 
g'  =  generation  coefficient  (cm”^sec“^) 

Equation  (3)  is  consistent  with  Russian  models  of  the  generation  rate 
(8:1529).  The  generation  coefficient,  g' ,  is  the  parameter  that  is 
usually  referenced  instead  of  the  generation  rate,  G' . 

The  second  rate  equation  is  the  net  ionization  of  traps,  or  more 
simply,  the  net  trapping.  For  master  model  one,  this  equation  models 
the  rate  of  trapping,  filling  of  traps,  and  deionization  of  traps.  The 
corresponding  equation  for  the  clipped  master  model  two  does  not  model 
the  filling  of  traps  or  the  deionization  of  traps.  Equations  (4)  and 
(5)  give  the  net  ionization  of  traps  rate  equations  for  master  models 
one  and  two,  respectively: 


dn.  ' 
'dT^ 


(4) 


dn.  ' 
~dF^ 


(5) 
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where 


n_'  =  density  of  ionized  traps  (cm“^) 

Nt'  =  density  of  traps  (cm~^) 

5"^'  =  trap  recombination  coefficient  (sec/cra^) 

Tj'  =  average  ionized  trap  lifetime  (sec) 

rj'  =  average  electron  drift  time  before  trapping  (sec) 

Equation  (5)  for  the  second,  or  clipped,  model  is  the  more  common  form 
in  most  Russian  publications  (8:1528,  12:879). 

The  final  rate  equation  models  the  time  rate  of  change  of  free 
electrons  in  the  conduction  band  of  the  crystal.  This  equation  is 
derived  from  electromagnetic  principles  which,  although  basic,  are  worth 
repeating  here  to  highlight  approximations.  We  start  with  two  of 
.Maxwell's  equations  given  in  Eqs  (6)  and  (7)  below: 

V  •  D'  -  V  •  eE '  -  p  '  ^ ^ ^ 

(7) 

The  "t",  "c" ,  "d" ,  and  "diff"  subscripts  on  the  current  density  terms, 
J',  refer  to  total,  conduction,  displacement,  and  diffusion  currents, 
respectively.  Ac  this  point,  approximations  are  made  for  the  total 
current.  In  some  papers,  Bryksin  and  others  ignore  everything  except 
conduction  current  (9:80).  In  most  works,  however,  only  the  diffusion 
current  is  neglected.  We  will  neglect  diffusion  current  in  both  our 
rna.srer  models  while  keeping  both  conduction  and  displacement  current. 
Bryksin  addresses  the  validity  of  this  approximation  in  one  of  his  1984 
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papers  (14;1491).  Under  most  conditions,  the  effects  of  diffusion 
current  are  negligible. 

The  development  of  the  final  rate  equation  continues  in  Eqs  (8)  to 
(11)  below.  Equation  (8)  is  simply  the  del  dot  product  of  both  sides  of 
Eq  (7),  with  total  current  approximated  by  conduction  current  and 
displacement  current: 


V  •  (oT/  + 


(9) 


D1 


V  •  J 


c 


0 


(10) 


-V  •  J 


c 


(11) 


Equation  (11)  has  used  Eq  (6)  and  the  fact  that  the  divergence  in  our 
one-dimensional  models  reduces  to  a  partial  derivative  with  respect  to 
the  longitudinal  dimension  of  the  crystal.  The  last  two  relationships 
needed  to  complete  the  development  are  given  in  Eqs  (12)  and  (13)  on  the 
following  page : 
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(12) 


Jc'  -  E' 


p  '  -  e(n,'  -  n,'  -  J]^') 


(13) 


where 

e  =  elementary  charge  (1.6022*10”^®  C) 

/j  =  mobility  of  electrons  (cm^/V»sec) 

Substituting  Eqs  (12)  and  (13)  into  Eq  (11)  and  rearranging  we  arrive  at 
the  fine"!  form  for  the  thi^'d  rate  equation.  Since  we  will  only  be 
dealing  in  one  dimen'  ,  the  vector  notation  will  be  omitted  for  the 
rest  of  the  repoi .  Equation  (14)  below,  therefore,  models  the  time 
rate  of  chap'^u  of  free  electrons  in  the  conduction  band: 


8n.  ' 

at' 


dn_ ' 
'dt^  * 


d 

dz' 


(pn^'E') 


(14) 


Equation  (14)  is  valid  for  both  master  models.  When  the  appropriate 
equations  for  the  first  two  terms  on  the  right  hand  side  of  Eq  (14)  are 
substituted,  however,  the  resulting  equations  for  the  two  master  models 
become  different. 

To  summarize,  the  three  rate  equations  for  the  comprehensive  master 
model  one  are  Eqs  (1),  (4-),  and  (14).  The  corresponding  rate  equations 
for  the  clipped  master  model  two  are  (2),  (5),  and  (14). 
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Normalization 


Equations  given  in  the  previous  section  were  dimensional.  To 
facilitate  further  mathematical  and  computer  calculations,  all  variables 
were  undimensionalized ,  or  normalized.  The  normalization  we  used  is 
based  on  the  initial  transient  time,  to',  which  is  derived  from  the 
crystal  thickness,  do',  the  electron  mobility  in  the  crystal,  /i,  and  the 
voltage  applied  to  the  transparent  electrodes,  Vq' .  In  this  report,  the 
initial  transient  time  will  also  be  referred  to  as  the  time  constant. 

The  initial  transient  time  is  simply  the  time  required  for  a  free 
electron  to  drift  from  the  cathode  to  the  anode  through  the  crystal 
under  the  influence  of  the  applied  electric  field,  assuming  no  trapping 
and  no  recombination  with  holes.  This  normalization  is  fairly  standard, 
however,  others  are  used.  For  example,  Bryksin  normalizes  to  "the 
characteristic  time  for  electron  trapping  by  a  trap"  in  one  of  his  1983 
papers  (7:686).  The  "characteristic  time  for  electron  trapping" 
corresponds  to  in  Eq  (5)  on  page  22.  The  normalization  most 
frequently  used  in  Russian  papers,  including  those  by  Bryksin,  is 
identical  to  the  normalization  presented  in  this  paper.  Table  1  on  the 
next  page  details  the  normalization  factors  and  subsequent  normalized 
variables.  Dimensional  parameters  such  as  e .  h,  e,  and  that 
represent  well  known  physical  constants  or  material  parameters,  do  not 
carry  the  dimensional  (')  notation  despite  the  fact  that  they  are 
dimensional.  Furthermore,  some  variables  such  as  I'  and  are  never 
explicitly  normalized  individually. 
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Table  1 


Normalization  Factors  and  Normalized  Variables 


NORMALIZATION  FACTORS 
do' 

to'  =  do'V(MVo') 

no'  =  €Vo'/(edo'2.102)* 

r 

I' 

Vo' 

M 

€  =  epfo 

^  r 

to  =  8.8542*10'12 

e  =  1.6022*10-i® 


DESCRIPTION 

Device  thickness  (cm) 

Initial  transient  time  (sec) 
Concentration  normalization  factor  (cm”^) 
Total  exposure  energy  (J/cm^) 

Write  light  intensity  (W/cm^) 

Applied  external  voltage  (Volts) 
Mobility  of  electrons  (cmVV'sec) 
Permittivity  (F/m) 

Relative  Permittivity 
Permittivity  of  free  space  (F/m) 
Elementary  charge  (C) 


Factor  of  10^  in  no'  is  a  conversion  factor  to  account  for  the 
difference  in  length  units  between  c  and  do'  . 


NORMALIZED  VARIABLES 

n+  =  n+'/no' 

n.  =  n_'/no' 

n^  =  ng  /no 

tis  ~  /^o 

Nd  =  Nd'/no' 

Nt  =  N,'/no' 

t  =  t'/to' 

^exp  ~  ^  /  ( 1  ^0  ) 

’’t  =  ’’t'/to' 

’’d  ’’d  /to 

z  =  z'/do' 
p  =  p' /(eno'  ) 

Q  =  a 'do' 

E  =  E'/(Vo'/do') 

J  =  J'/(eno'do'/to' ) 
Cd  =  rd'/(no'to'  ) 

Tt  “  Tt  /(r>o  to  ) 

G  =  G'to' 

g  =  g'/(no'do'/to'  ) 


DESCRIPTION 

Density  of  ionized  donors 

Density  of  ionized,  or  filled,  traps 

Density  of  free  electrons 

Injection  coefficient  of  cathode 

Density  of  donors 

Density  of  traps 

Time 

Exposure  time 

Average  ionized  trap  lifetime 
Average  electron  drift  time 
Longitudinal  distance  or  crystal  depth 
Charge  density 
Absorption  coefficient 
Electric  field  intensity 
Current  density  (J^,  Jg,  J^,  or  Jdiff) 
Donor  recombination  coefficient 
Trap  recombination  coefficient 
Generation  rate 
Generation  coefficient 
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After  normalization,  the  three  rate  equations  for  master  model  one 


become 


dn^ 

'Wt 


g-^  {N^  -  nj  exp(-az) 


(15) 


dn  {N^  -  n_)  n 

^=5—  ■  ^ -  -  — 

dc  ^  Ct  T 


(16) 


dc 


dn, 

~dF 


dn, 

dc 


(17) 


Similarly,  the  normalized  rate  equations  for  the  clipped  master 
model  two  are 


dn, 

'dF 


g-^  (N^  -  n.)exp(-az) 


(18) 


dn_  n^ 

dc 


(19) 


dng  dn, 

dc  dc 


dn. 

~dC 


(17) 
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Computer  Modelling 


Before  the  rate  equations  can  be  solved  numerically  with  a  finite 
difference  scheme ,  they  must  be  converted  to  a  discrete  form  and  initial 
and  boundary  conditions  must  be  specified.  To  accomplish  this,  the 
crystal  is  divided  longitudinally  into  100  discrete  and  equal  bins. 
Accuracy  will  improve,  of  course,  with  an  increased  number  of  bins.  Our 
differential  equation  solver  software,  Mathlib  DEQSOLVE  by  Innosoft 
International,  Inc.,  would  accept  up  to  200  bins.  The  Russians 
apparently  used  201  bin.s  (8:1530). 

The  time  step  over  which  a  solution  is  attempted  must  also  be 
specified.  The  magnitude  of  the  time  step  was  manipulated  by  specifying 
an  absolute  tolerance  value  to  the  differential  equation  solver 
software.  The  software  then  internally  determined  the  appropriate  time 
step  to  ensure  the  absolute  tolerance  value  was  not  exceeded.  Early  on, 
we  determined  that  an  absolute  tolerance  value  of  10~^  resulted  in 
excessive  computation  times  and  did  not  significantly  improve  accuracy. 
As  a  result,  a  tolerance  value  of  10"^  was  used  for  most  runs.  The 
Russians  use  a  finite  difference  scheme  with  FORTRAN  EC.  The  specifics 
of  the  program  are  not  provided,  however,  they  profess  an  accuracy 
better  than  10'^  in  previous  reports  (8:1530). 

Since  the  discrete  nature  of  the  time  step  was  made  transparent  co 
the  user  by  the  software,  the  first  two  rate  equations  are  inherently 
discrete  and  require  only  minor  modifications  in  order  to  be  acceptable 
to  the  software.  The  third  rate  equation,  however,  had  to  be  modified 
as  follows : 
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AJ^e  _  ^  A(73g£’)  (20) 

At  At  AZ 


Using  secondary  subscripts  to  denote  the  bin  number  and  the  fact  that, 
after  normalization,  1/Az  equals  nbin,  the  last  term  in  Eq  (20)  is  more 
conveniently  expressed  by 


A(n^g) 

AZ 


(21) 


After  the  rate  equations  have  been  made  discrete  and  appropriate 
parameters  entered,  the  solution  can  proceed.  The  initial  values  for 
density  of  ionized  donors,  ionized  traps,  and  free  electrons  are  set  to 
zero,  and  cathode  and  anode  boundary  conditions,  appropriate  to  the  type 
of  FRIZ,  are  entered.  For  the  DI  and  SI  FRIZ,  the  boundary  value  for 
neo  is  set  to  zero  since  no  injection  is  allowed  through  the  cathode. 

For  the  conducting  FRIZ,  the  value  of  ngo  must  reflect  injection.  To 
accomplish  this,  we  have  used  the  same  model  as  Bliznetsov  reports  in 
the  1988  paper  (6:375).  The  dimensional  form  for  this  equation  is 
actually  specified  by  Bliznetsov,  however,  normalization  does  not  change 
the  equation's  form  as  can  be  seen  by  Eq  (22)  on  the  next  page: 
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n^g{  t)  -  n 


1 


gpto) 

EoU) 


(22) 


where 

neo(t)  =  electron  density  at  the  cathode  at  time  t 
ng  =  cathode  injection  coefficient 
Eq(0)  =  electric  field  at  cathode  at  time  0 
EgCt)  =  electric  field  at  cathode  at  time  t 

The  value  of  at  the  anode,  or  ng^oQ.  is  automatically  calculated 
by  the  routine.  For  a  DI  PRIZ,  however,  the  term  neiooEioo(^)  =  -^cioo 
is  set  to  zero  since  there  can  be  no  current  through  the  anode.  The 
various  values  for  Eo(t),  the  cathode  field,  and  Ej(t),  where  the 
subscript  j  refers  to  the  bin  number,  are  calculated  with  Eqs  (23)  and 
(24)  below: 


n^in 


(23) 


(24) 


Equations  (23)  and  (24)  are  derived  from  manipulations  of  the  discrete 
form  of  Eq  (6).  The  derivation  of  Eqs  (23)  and  (24)  forces  the  integral 
of  the  electric  field  across  the  length  of  the  device  to  equal  -Vq  for 
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all  time.  The  validity  of  the  derivation  is  tested  during  all 
calculations  by  simultaneously  calculating  the  net  accumulated  charge 
density  within  the  crystal  as  time  progresses.  A  DI  PRIZ  will  always 
have  a  zero  net  accumulated  charge  density,  since  no  charges  can  enter 
or  leave  the  crystal,  however,  the  SI  and  conducting  PRIZ  will  develop  a 
net  charge  over  time. 

The  only  parameter  that  remains  to  be  specified  is  the  exposure 
time.  The  exposure  time  must  be  entered  in  normalized  time  units 
according  to  Eq  (25) : 


t 


exp 


(25) 


where 

^exp  =  normalized  exposure  time 
=  total  exposure  energy  (J/cm^) 

I'  =  write  light  intensity  (Watts/cm^) 

Given  appropriate  entries  for  all  the  parameters  specified  in  the 
paragraphs  above,  the  software  will  compute  ionized  donor  density, 
ionized  trap  density,  and  free  electron  density  versus  time  and 
distance.  Using  this  information,  the  net  charge  density  can  be 
computed,  as  well  as  the  internal  electric  field  and  current  density, 
versus  position  and  time.  The  computation  time  for  one  run  is  anywhere 
from  15  minutes  to  several  hours  depending  on  exposure  time,  write  light 
intensity,  and  numerous  other  factors.  For  example,  since  the  number  of 
differential  equations  that  must  be  solved  for  each  time  step  is  nbin 
times  three,  the  computation  times  increase  with  an  increase  in  nbin. 
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IV .  Detailed  Analysis  of  Soviet  Journal  Articles 


Discussion  and  Interpretation 

The  two  Soviet  papers  we  specifically  address  in  this  thesis  are 
"Investigation  of  Reciprocity  in  the  PRIZ  Space-Time  Light  Modulator", 
by  Bliznetsov,  Kuz'min,  and  Khomenko  (1988  paper)  and  "Injection 
Mechanism  of  Dynamic  Image  Selection  in  PRIZ  Space-Time  Light 
Modulators",  by  Bliznetsov,  Bryksin,  Korovin,  Miridonov,  and  Khomenko 
(1987  paper) .  Most  of  our  information  on  Soviet  analytical  modelling 
techniques  is  derived  from  previous  works  by  Bryksin.  Bryksin  is 
referenced  in  the  1987  and  1988  papers,  as  are  Petrov,  Shlyagin, 
Astratov,  and  others.  Bliznetsov  does  not  appear  as  an  author  in  any  of 
the  additional  references  we  have  used.  His  work  on  the  conducting  PRIZ 
appears  to  be  limited  to  the  1987  and  1988  papers,  and  the  extent  to 
which  he  has  chosen  unique  nomenclature  or  modelling  techniques  is 
unclear.  Both  papers  contain  data  that  is  not  fully  specified  and 
subject  to  some  interpretation.  The  experimental  data  presented  in  both 
papers  is  especially  nebulous  and,  as  a  result,  we  experienced  greater 
success  in  modelling  the  Russian's  analytical  data  than  their 
experimental  data.  Even  the  Russians,  however,  claim  to  analytically 
model  only  the  "qualitative"  aspects  of  their  own  experimental  data 
(6:374) . 

In  the  1988  paper,  the  authors  initially  present  experimental  data 
for  the  operation  of  a  conducting  PRIZ  over  the  range  of  writing  light 
intensities  from  30  /iW/cm^  to  300  mW/cm^.  Their  PRIZ  consists  of  a 
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0.45  mm  thick  Bii2Si02o  (BSO)  crystal,  with  electrodes  biased  by  a  two 
kilovolt  DC  power  supply.  The  image,  a  20  fim  wide  strip,  is  written 
with  a  helium-cadmium  laser  (A  =  441  run).  The  imaged  strip  produces 
transverse  electric  fields  within  the  PRIZ,  and  a  subsequent  response. 
The  image  is  read  with  a  helium-neon  laser  (A  -  633  nm)  with  an 
intensity  of  10  /iW/cm^. 

As  a  basis  for  our  interpretation,  we  assumed  that  the  response 
measured,  and  experimentally  reported,  by  the  Russians  was  the 
instantaneous  output,  or  throughput,  of  the  read  light  for  a  given 
condition  of  exposure.  In  other  words,  no  time  delay  was  used  to 
measure  the  response  of  the  PRIZ  to  an  applied  stimulus.  The  impact  of 
time  delay  on  response  will  be  covered  in  detail  in  subsequent  sections. 
At  this  point,  we  will  simply  mention  that  the  finite  time  required  for 
photogenerated  electrons  to  drift  under  the  influence  of  the  electric 
field  within  the  crystal  results  in  a  delay  of  response  to  a  given 
stimulus.  Data  can  be  collected  by  either  allowing  for  or  ignoring  this 
delay.  We  believe  rather  confidently  that  the  Russians  have  chosen  the 
latter  option.  Furthermore,  at  some  point,  we  must  link  the  results 
from  our  one  dimensional  analytical  model  to  the  response  data  from  the 
multidimensional  experiments.  The  Russians  never  attempted  this,  and 
instead,  relied  on  qualitative  explanations  to  link  their  experimental 
and  analytical  data.  Using  these  qualitative  explanations,  we  attempted 
to  quantify  the  response  produced  by  our  various  one  dimensional  charge 
distributions.  This  consisted  of  predicting  the  response  produced  when 
a  one-dimensional  charge  distribution  is  placed  adjacent  to  a  net  zero 
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charge  distribution,  resulting  in  transverse  electric  fields  that 
modulate  the  read  light.  The  criterion  we  chose  to  quantify  the 
transverse  fields  is  crude,  yet  consistent  with  qualitative 
explanations.  We  will  elaborate  on  the  details  of  our  criterion  later 
in  this  section.  With  these  overriding  assumptions  in  mind,  we  continue 
with  our  detailed  analysis  of  the  Russian  papers. 

First,  the  authors  experimentally  establish  a  region  over  which  the 
principle  of  reciprocity  holds.  The  principle  of  reciprocity  asserts 
that  the  response  of  the  FRIZ  is  independent  of  write  light  intensity 
for  a  given  exposure,  where  the  exposure  is  given  by  the  product  of  the 
write  light  intensity  and  time  of  illumination.  They  define  response 
somev;hat  ambiguously  as  "the  total  intensity  3?  of  the  readout  image  at 
various  reading-light  intensities."  (6:374)  The  term  "reading-light 
intensities"  is  certainly  intended  to  be  writing-light  intensities, 
however,  further  interpretation  of  their  definition  is  difficult. 
Equatl''.,  (26)  expresses  the  principle  of  reciprocity,  as  reported  in  the 
1988  paper,  for  values  of  exposure  less  than  20  /uJ/cm^: 

(26) 


where 


S  =  FRIZ  response 

=  exposure  (J/cm^) 

I'  =  write  light  intensity  (W/cm^) 
^exp"  “  exposure  time  (sec) 
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Figure  3,  below,  approximates  experimental  data,  presented  in 
Figure  One  of  the  1988  paper,  backing  up  the  authors'  claims  concerning 
reciprocity.  Specifically,  we  can  see  from  Figure  3  that  there  is  a 
region  where  the  different  intensity  curves  coincide,  indicating 
reciprocity.  We  can  also  see  that  outside  this  reciprocity  region,  for 
a  given  exposure,  the  higher  intensity  curve  produces  a  higher  peak 
response.  Furthermore,  we  can  see  that  the  peak  response  (Umax)  occurs 
at  higher  exposures  (^max)  higher  intensities.  Unfortunately,  the 

original  graph  (Figure  One,  1988  Paper)  is  devoid  of  any  numerical  data 
and  we  can  not  be  sure  what  intensities  have  been  drawn.  Despite  thi.s, 
the  authors'  assert,  "As  shown  schematically  in  Fig.  One,  at  low 
exposures  <  20  ^J/cm^)  the  functions  coincide  to  within  30%  for 

all  values  of  the  writing-light  intensity  in  the  range  from  3 • 10"^  to 
3 • 10'^  W/cm^."  (6:374)  "Fig.  One"  referred  to  above  is  our  Figure  3. 

The  authors  present  a  wealth  of  information  on  the  trends  in 
^max-  This  information  is  reproduced  in  Figure  4. 

Having  provided  experimental  data  on  reciprocity,  the  authors 
proceed  to  explain  why  reciprocity  holds  below  a  certain  value 
(<  20  ^J/cm^)  and  fails  above  this  value.  In  developing  their 
explanation,  they  assert  that  the  response  of  the  PRIZ  is  due  to  the 
formation  of  a  positive  space  charge  near  the  cathode.  The  magnitude 
and  extent  of  this  space  charge  is  determined  by  two  competing  factors. 
The  first  factor  is  the  photoionization  of  donors  near  the  cathode, 
producing  electrons  that  drift  toward  the  anode  due  to  the  electric 
field.  This  effect  increases  the  positive  space  charge  near  the 
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S  iJJ/cm^ 


Figure  3.  Typical  Dependencies  of  the  Intensity  IR  of  the  Readout  Image 
on  the  Exposure  Energy  for  Various  Intensities  I'  of  the  Writing 
Light:  I'(l)  «  I' (2)  «  I' (3).  Taken  from  Figure  1,  88  Paper  (6:374) 


Figure  4.  Experimental  Dependencies  of  the  Maximum  Intensity  R^ax 
Readout  Image  and  the  Corresponding  Exposure  Energy  the  Writing- 

Light  Intensity  I'.  Taken  from  Figure  2,  88  Paper  (6:375) 


X  arb.  units 


cathode.  The  second  factor  is  the  injection  of  electrons  through  the 
cathode.  As  previously  discussed,  the  injection  of  electrons  will 
increase  as  the  cathode  electric  field  increases  (becomes  more 
negative).  The  injection  of  electrons  tends  to  reduce  the  positive 
space  charge  near  the  cathode.  The  authors  use  the  interaction  of  these 
two  competing  factors  to  explain  both  the  reciprocal  and  nonreciprocal 
regions . 

To  explain  reciprocity,  the  authors  state,  "In  the  initial  stage  of 
the  writing,  when  the  charge  density  is  low,  the  injection  current  is 
low  in  comparison  with  the  photocurrent  and  does  not  have  any 
appreciable  effect  on  the  charge  accumulation."  (6:374)  The  sentences 
in  the  1988  paper  preceding  the  one  above  imply  that  this  statement  only 
applies  under  high  intensity  conditions.  For  high  intensities  and  low 
exposures,  the  total  exposure  time  corresponds  to  the  "initial  stage" 
mentioned  above.  Injection,  therefore,  is  not  as  great  a  factor  at  high 
intensities,  and  low  exposures,  as  it  is  for  low  intensities.  As  we 
will  soon  see,  for  intensities  of  500mW/cm^,  and  exposure  values  within 
the  reciprocity  region,  injection  is  virtually  nonexistent.  This  occurs 
because  the  characteristic  exposure  times,  t^^p,  at  high  intensities  and 
low  exposures,  are  less  than  the  drift  time  of  the  photoelectrons.  As  a 
result,  the  space  charge  within  the  crystal  can  not  fully  develop  and 
the  cathode  field  remains  near  its  initial  value,  allowing  little 
injection.  At  low  intensity  and  low  exposure,  the  characteristic 
exposure  times  are  long  enough  to  allow  significant  injection. 
Presumably,  the  relative  lack  of  injection  at  high  intensities  and  low 
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exposures  explains  why  the  response  of  the  PRIZ  is  reciprocal  at 
exposures  less  than  20  /iJ/cm^,  The  preceding  explanation  is 
unsatisfying,  as  it  intuitively  appears  to  argue  against,  rather  than 
for,  reciprocity  at  low  exposures.  In  other  words,  at  low  intensity  and 
low  exposure,  the  buildup  of  positive  space  charge  near  the  cathode  is 
counteracted  to  some  degree  by  injection.  At  high  intensity  and  the 
same  low  exposure,  however,  the  buildup  of  positive  space  charge  near 
the  cathode  proceeds  unimpeded  by  injection.  If  response  is  indeed  due 
to  the  density  of  positive  space  charge  near  the  cathode,  as  the  authors 
assert,  then  the  higher  intensity  light  should  result  in  a  greater  net 
positive  space  charge  and,  therefore,  a  higher  response. 

The  authors'  explanation  for  the  failure  of  reciprocity  at  higher 
exposures  is  a  little  more  satisfying,  even  though  it  seems  to  be 
inconsistent  with  their  previous  explanation  for  the  existence  of 
reciprocity  at  low  exposures.  In  their  words,  "the  main  factor 
determining  the  deviation  from  the  reciprocity  law  for  the  PRIZ 
modulator  is  the  injection  of  electrons  through  the  contact,  which 
causes  the  photoinduced  positive  charge  to  be  partially  compensated." 
(6:376).  As  we  have  already  seen,  they  use  a  very  similar  argument  to 
explain  the  reciprocity  region  at  low  exposures.  Another  statement  the 
authors  make,  in  explaining  the  failure  of  reciprocity,  serves  as  the 
basis  for  our  attempts  to  explain,  in  one  dimension,  an  essentially  two- 
dimensional  process.  Again,  in  their  words,  "Since  an  increase  in  the 
density  of  positive  charge  brings  about  an  increase  in  the  amplitude  of 
modulation  of  the  reading  light,  the  maximum  intensity  of  *^^6 
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readout  image  increases  at  high  intensities  of  the  writing  light." 
(6:374)  Here,  the  authors  have  connected  their  unspecified  response  (R) 
to  a  parameter  we  can  easily  model  in  one  dimension,  the  positive  charge 
density.  This  connection  forms  the  basis  for  our  efforts  to  unravel  the 
apparent  inconsistencies  in  the  1988  paper  concerning  reciprocity. 

Next,  the  1988  paper  mentions  a  mathematical  model  for  the  SI  and 
conducting  PRIZ  that  the  authors  assert  will  verify  the  "qualitative 
explanations"  given  for  reciprocity  (6:374).  The  model  is  said  to  be 
developed  from  three  previous  works  by  Bryksin,  which  we  also  used  in 
our  models.  Unfortunately,  the  three  Bryksin  references  were  not 
consistent.  One  of  the  references  did  not  model  trapping  at  all,  and 
another  neglected  displacement  current  in  the  formula  for  total  current, 
presumably  in  addition  to  neglecting  diffusion  current.  Furthermore, 
one  of  the  references  models  recombination  of  donors  and  detrapping 
while  the  others  do  not.  In  the  1988  paper's  model,  the  authors 
specifically  claim  to  take  into  account  the  ionization  and  depletion  of 
donors,  the  trapping  of  electrons,  the  decreasing  intensity  of  light 
versus  depth  of  penetration  (absorption  coefficient  modelling) ,  and  the 
injection  of  electrons  through  the  cathode  at  a  rate  dependent  upon  the 
electric  field  at  the  cathode.  The  authors  never  provide  the  equations 
for  their  entire  model,  although  they  do  specify  the  form  for  the 
injection  as  previously  specified  in  Eq  (22).  The  only  hint  they  give 
as  to  most  of  the  specific  material  parameters  they  use  is  that  the 
values  "corresponded  to  the  conditions  of  the  experiment."  (6:374) 
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The  Russians  first  use  their  model  to  compare  the  net  space  charge 
density  in  an  SI  and  conducting  PRIZ  under  low  intense  ty  illumination 
conditions  (I'  =  50  /iW/cm^)  .  They  plot  data  for  net  space  charge 
density,  p,  versus  longitudinal  distance,  z' ,  for  exposures  of  0.4,  2, 
and  4  pj/cm^.  Interestingly,  they  do  not  report  any  data  for  exposures 
outside  the  established  reciprocity  region.  This  would  seem  to  be  a 
logical  step  if  they  are  trying  to  analytically  establish  injection  as 
the  cause  for  the  failure  of  reciprocity  at  higher  exposures.  Even  more 
mysteriously,  some  of  their  data  appears  to  show  that  reciprocity 
analytically  fails  where  it  experimentally  succeeded.  Figures  5  and  6 
are  normalized  reproductions  of  their  plots  for  the  SI  and  conducting 
PRIZ,  respectively.  Though  the  graphs  appear  similar,  a  quick  glimpse 
at  the  scales  will  show  that  the  conducting  PRIZ  produces  net  positive 
space  charges  that  are  less  than  half  as  high  and  wide  as  the 
corresponding  values  for  the  SI  PRIZ.  As  predicted  earlier,  at  low 
intensity  and  low  exposure,  injection  reduces  the  positive  space  charge 
region.  As  they  state,  "at  low  writing-light  intensities  injection  of 
electrons  from  the  contact  appreciably  reduces  both  the  density  of  the 
positive  charge  and  the  thickness  of  the  layer  that  it  occupies  in  the 
crystal."  They  further  state,  "Injection  under  these  conditions  thus 
decreases  the  magnitude  of  the  STLM  response."  (6:375)  STEM  is  an 
acronym  for  space-time  light  modulator.  With  Figures  5  and  6,  the 
authors  have  established  that  injection  is  significant  at  low 
intensities,  for  exposures  in  the  reciprocity  region.  If  they  can 
subsequently  show  that  injection  is  insignificant  at  high  intensities, 
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Figure  5.  Net 
I'  =  50  ^(W/cm^, 
3a,  88  Paper  (6: 


Figure  6.  Net  S' 
with  Hg  =  0.1,  I 
Taken  from  Fieur 


for  exposures  within  the  reciprocity  region,  then  some  of  their  earlier 
claims  are  justified  even  though  reciprocity  is  still  not  clearly 
established.  To  fully  demonstrate  reciprocity,  we  must  plot  net  space 
charge  densities  for  the  conducting  PRIZ  at  various  intensities  and 
exposures,  within  and  without  the  region  of  reciprocity.  A  comparison 
of  conducting  PRIZ  reciprocity  data  with  SI  PRIZ  reciprocity  data  should 
also  serve  to  highlight  the  specific  effects  of  injection  on 
reciprocity. 

Next,  the  Russians  present  data  for  a  conducting  PRIZ  under  high 
write  light  intensity  (I'  =  500  mW/cm^)  for  exposures  of  4,  8,  12,  and 
24  pj/cm^.  Again,  a  seemingly  strange  choice  has  been  made,  as  the 
intensities  are  well  outside  the  experimental  data,  and  the  exposures 
are  all  less  than  or  comparable  to  the  reciprocity  region  boundary.  The 
Soviet  data  is  presented  in  normalized  form  in  Figure  7.  The  authors 
assert  that,  "As  the  intensity  of  the  writing  light  increases,  the  range 
of  exposures  at  which  injection  has  only  insignificant  effect  becomes 
broader."  (6:375)  Though  the  data  is  absent,  the  authors  reinforce  the 
assertion  above  by  stating  that  the  SI  PRIZ  data  for  500  mW/cm^  differs 
by  less  than  one  percent.  If  this  is  true,  then  the  basis  for  some  of 
the  authors'  claims,  as  described  in  the  previous  paragraphs,  is 
established.  Still,  however,  no  direct  comparison  is  done  for  space 
charge  densities  at  various  intensities  for  specific  exposure  values. 

One  comparison  is  possible  by  plotting  the  4  pj/cm^  exposure  curves  from 
Figures  6  and  7  on  the  same  chart.  The  result  of  this  exercise  is  shown 
in  Figure  8.  The  two  space  charge  densities  in  this  plot  are 
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Figure  7.  Net  Space  Charge  Density,  p,  versus  t  for  a  Conducting  PRIZ 
with  ns  =  0.1,  I'  =  500  mW/cm^,  r"  =  4  (1) ,  8  (2)  ,  12  (3),  and  24  pJ/cm^ 
(4).  Taken  from  Figure  4,  88  Paper  (6:375) 
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Figure  8.  Comparison  of  Net  Space  Charge  Densities,  p,  for  ng  -  0.1, 

=  4  pJ/cm^,  I'  =50  pW/cm^  (1)  and  I'  -  500  mW/cm^  (2).  Taken  from 
Figures  3b  c^nd  4,  88  Paper  (6:375) 


significantly  different,  even  though  all  experimental  data  and 
explanations  given  up  until  now  indicate  that  they  should  be  similar,  or 
at  least  produce  similar  transverse  fields,  indicating  reciprocity.  The 
unexpected  differences  in  the  two  space  charge  distributions  prompted  us 
to  develop  a  simple  criterion  by  which  to  gauge  the  transverse  fields, 
and  subsequent  response,  produced  by  a  given  one-dimensional  space 
charge  distribution.  This  criterion  is  detailed  in  the  next  paragraph. 
An  even  more  puzzling  aspect  highlighted  by  Figure  8  is  the  fact  that 
the  high  intensity  illumination  produces  a  smaller  and  narrower  space 
charge  region  than  the  low  intensity  illumination.  This  is  opposite 
from  the  effect  one  would  expect  from  previous  explanations. 

Since  the  authors  never  explain  the  exact  connection  between  net 
positive  charge  distribution  and  response,  we  must  establish  some 
definite  criterion  for  determining  whether  two  one-dimensional  space 
charge  plots  produce  the  same  response,  within  a  30  percent  tolerance  as 
specified  in  the  1988  paper  (6:37A).  Remembering  that  FRIZ  response  is 
due  to  transverse  electric  fields  within  the  crystal,  we  need  a 
criterion  that  consistently  represents  the  response  produced  by  the 
transverse  electric  fields  that  result  from  placing  our  one-dimensional 
space  charge  distribution  adjacent  to  a  region  with  a  uniformly  zero 
space  charge  distribution.  The  boundary  between  our  space  charge 
distribution  and  the  adjacent  net  zero  distribution  corresponds  to  the 
edge  of  the  slit  that  is  imaged  onto  the  device  by  the  write  light. 

Since  the  authors  assert  that  the  response  is  due  to  both  the  magnitude 
of  the  density  of  positive  charge  and  the  thickness  of  the  positive 
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charge  layer,  the  criterion  we  will  use  is  the  area  under  the  positive 
portion  of  the  net  space  charge  curve  (6:375).  If  the  areas  under  two 
different  curves  are  within  30  percent  of  each  other,  we  will  assume 
that  they  produce  responses  within  30  percent  of  each  other.  Using  this 
criterion  alone,  the  curves  in  Figure  8  do  not  represent  reciprocity. 

The  crudeness  of  our  criterion  does  not  allow  a  total  rejection  of 
reciprocity  in  this  case,  or  in  many  other  cases  which  we  will  analyze. 
Our  results  in  the  next  section  are  just  as  much  a  test  of  the  validity 
of  our  criterion  as  they  are  a  study  of  reciprocity.  Nevertheless,  we 
felt  some  value  could  be  obtained  from  even  a  crude  attempt  to 
characterize  the  response  produced  by  a  given  one-dimensional  space 
charge  distribution.  Perhaps  the  most  obvious  weakness  of  our  criterion 
is  the  total  disregard  of  the  transverse  fields  produced  by  the  region 
of  net  negative  space  charge  density.  Other  weaknesses  will  be 
discussed  in  the  next  section. 

In  concluding  the  paper,  the  authors  assert  that  the  response  of 
the  FRIZ  lags  the  illumination  of  the  writing  light  for  conditions  of 
high  intensity  illumination  and  corresponding  low  exposure  times.  This 
will  occur  any  time  the  exposure  time  is  of  the  same  magnitude  or  less 
than  the  drift  time  of  an  electron  that  has  been  photogenerated.  This 
final  comment  sheds  a  little  more  light  on  the  authors'  claims 
concerning  reciprocity.  Accepting  Figure  8  as  the  true  depiction  of  the 
space  charge  distributions  for  high  and  low  intensity  illumination,  in 
spite  of  previously  mentioned  discrepancies,  the  effects  of  lag  time  can 
be  used  to  account  for  the  apparent  lack  of  reciprocity.  Lag  time 
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explains  that  the  space  charge  distribution  for  the  high  intensity 
illumination  is  not  as  great  since  charge  separation  has  only  begun. 
Perhaps,  if  the  response  is  measured  after  the  FRIZ  has  had  time  to 
develop,  the  high  intensity  space  charge  curve  should  shift  outward  and 
indicate  reciprocity.  The  shift,  or  decay,  of  space  charge 
distributions  following  turn-off  of  the  write  light  will  be  investigated 
in  the  next  section.  The  authors  appear  to  indicate  that  the  effects 
of  lag  time  and  lack  of  injection  at  high  intensity  exposures  counteract 
each  other  and  produce  a  combined  effect  that  demonstrates  reciprocity, 
relative  to  the  corresponding  low  intensity  situation. 

The  1987  paper  ultimately  investigates  two-dimensional  aspects  of 
dynamic  image  selection  (DIS) ,  an  operating  mode  of  the  PRIZ  in  which 
stationary  inputs  are  suppressed  while  moving  inputs  are  imaged.  As 
they  state  in  the  preamble,  "It  is  shown  that  the  dominant  role  in  the 
onset  of  the  response  of  the  modulator  to  the  turn-off  of  the  writing 
light  is  played  by  a  mechanism  involving  the  spatially  nonuniform 
injection  of  electrons  into  the  volume  of  the  crystal  through  the 
negative  electrode-crystal  contact."  (5:750)  In  arriving  at  this 
conclusion,  the  authors  present  experimental  data  on  the  time  evolution 
of  diffraction  efficiency  for  writing  light  intensities  of  30  nV/cm^ 
and  for  various  values  of  the  turn-off  time,  or  exposure.  One  of  the 
primary  features  of  this  data  is  the  increasing  of  diffraction 
efficiency  with  exposure  for  low  exposures.  This  is  consistent  with 
data  in  the  1988  paper  if  one  accepts  diffraction  efficiency  and 
response  as  similar  parameters.  The  authors  then  assert  that  for  higher 
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exposures,  the  curves  immediately  plunge  to  zero  following  turn-off  of 
the  write  beam  and  then  proceed  to  regenerate  themselves.  Depending 
upon  interpretation,  this  data  modifies  or  perhaps  even  contradicts  data 
in  the  1988  paper.  Unfortunately,  limitations  of  our  one— dimensional 
models  did  not  permit  adequate  examination  of  these  apparent 
contradictions  or  the  phenomena  of  regeneration. 

Specific  Modelling  Goals 

Our  modelling  goals  are  divided  into  four  major  categories,  which 
will  be  presented  separately  in  section  V.  These  categories  are 
reproduction  of  low  intensity  data  (Figures  5  and  6) ,  reproduction  of 
high  intensity  data  (Figure  7),  demonstrations  of  reciprocity  (Figures  3 
and  8) ,  and  demonstrations  of  self  erasure  (Figure  3) . 

Prior  to  the  first  category,  reproduction  of  low  intensity  data,  we 
will  include  a  validation  of  our  two  master  models.  We  will  accomplish 
this  using  the  intensity  and  exposure  data  specified  for  Figures 
5  and  6.  In  addition  to  the  data  presented  in  these  figures,  however, 
we  will  also  present  data  for  the  DI  PRIZ,  and  qualitatively  analyze  the 
DI  data's  consistency  relative  to  the  SI  and  conducting  PRIZ  data. 

After  validating  the  master  models,  we  will  attempt  to  reproduce 
Figures  5  and  6. 

In  the  second  category,  reproduction  of  high  intensity  data,  we 
will  attempt  to  reproduce  Figure  7.  We  will  also  verify  that  the 
corresponding  plots  for  an  SI  PRIZ  differ  by  less  than  one  percent  as 
the  authors  assert.  Attempts  will  also  be  made  to  verify  the  statement 
that,  "as  the  intensity  of  the  writing  light  increases,  the  range  of 
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exposures  at  which  injection  has  only  insignificant  effect  becomes 
broader."  (6:375) 

In  the  third  category,  demonstrations  of  reciprocity,  we  will 
present  a  full  range  of  data,  similar  to  that  contained  in  Figure  8,  to 
validate  the  experimentally  established  reciprocity  region.  This  will 
include  runs  at  a  wide  range  of  intensities  and  exposures  for  both  the 
SI  and  conducting  PRIZ.  This  should  shed  light  on  the  assertion  that 
the  response  of  the  PRIZ  is  due  to  the  net  positive  space  charge  density 
near  the  cathode . 

The  fourth  category  will  investigate  the  phenomena  of  self  erasure. 
Self  erasure  is  demonstrated  in  Figure  3  by  the  fact  that  response 
eventually  appears  to  decrease,  with  increased  exposure,  for  all 
intensities.  Once  again,  since  PRIZ  response  is  a  multi-dimensional 
phenomenon,  we  must  intuitively  speculate  on  the  nature  of  the 
transverse  fields  produced  when  our  one— dimensional  space  charge 
distribution  is  placed  adjacent  to  a  net  zero  space  charge  distribution. 
Using  this  analogy,  we  should  be  able  to  shed  some  light  on  the 
Russians'  assertions  concerning  the  correspondence  of  PRIZ  response  to 
the  net  positive  space  charge  density  near  the  cathode. 
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V .  Results 


Model  Validation 

Before  applying  our  two  master  models,  we  validated  them  by 
analyzing  the  data  they  produced  for  an  SI,  conducting,  and  DI  PRIZ 
under  conditions  corresponding  to  Figures  5  and  6  (I'  =  50  /iW/cm^, 

?“  =  4  ^J/cm^,  ng  =  0  for  SI  and  DI ,  =  0.1  for  conducting).  Three 

validation  runs  were  made  for  each  of  the  PRIZ  types.  These  runs 
corresponded  to  master  model  one  with  parameter  set  one,  master  model 
two  with  param.eter  set  one,  and  master  model  two  with  parameter  set  two. 
The  two  parameter  sets  are  listed  in  Table  2  with  appropriate  references 
to  justify  respective  values.  The  primary  differences  between  the  two 
parameter  sets  are  values  for  electron  mobility,  density  of  donors,  and 
absorption  coefficient.  Parameter  set  one  contains  values  that  appear 
to  predominate  in  most  of  the  Russian  references.  Many  of  the  values  in 
parameter  set  two  are  taken  from  references  by  Peltier,  Micheron, 
Horwitz,  and  Corbett  (30:3683,  25:353).  Two  parameters  that  require 
additional  explanation  are  the  donor  recombination  coefficient  and  the 
trap  recombination  coefficient.  The  value  for  the  trap  recombination 
coefficient,  is  obtained  from  experimental  data  for  times  such  that 

n_  «  Nj.  Thus,  its  value  is  estimated  by  multiplying  the  average 
electron  drift  time,  by  the  density  of  traps,  Nj .  In  effect,  this 

transformation  modifies  the  average  electron  drift  time  to  account  for 
the  depletion  of  traps  term  in  master  model  one.  The  slight  discrepancy 
in  the  value  for  ("t  reported  in  Table  2,  362  reported  versus  329  actual. 
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Table  2 


Parameter 

Sets 

Parameter  Set  1 

Parameter  Set  2 

Additional 

Parameter 

Value  (Reference) 

Value  (Reference) 

References 

nbin 

100 

100 

n 

1  (14;1A91) 

1  (38:1677) 

y’  (Hz) 

6.8-10^'* 

6.8-103'- 

(6:374) 

Vo'  (V) 

2000 

2000 

(6:374) 

do'  (cm) 

0.045 

0.045 

(6:374) 

jj.  (cm^/V*sec) 

0.01  (14:1487) 

0.03  (25:360) 

(8:1532) 

to'  (sec) 

1.0125*10-'* 

3.375-10‘3 

^  r 

50  (8:1532) 

56  (25:360) 

(14:1491) 

€  (F/m) 

4.427*10"i° 

4.958-10-30 

e  (C) 

1.6022«10'13 

1.6022-10-39 

no'  (cm~^) 

2.729-10^3 

3.056-1033 

q'  (cm’*) 

40  (38:1673) 

45  (11:1689) 

a 

1.8 

2.025 

Nd 

1-103  (8:1532) 

3.27-103  (30:3684) 

Nt 

366 

- 

(30:3684) 

’’d 

0.9 

0.9 

(11:1689) 

’’t 

1-103 

- 

(3:1585) 

Cd 

3.62-103 

- 

Tt 

3.62-102 

- 

ns 

0.1 

0.09 

(6:375) 

g  (I '=30  /iW/cm2) 

0.006 

0.0018 

g  (I '=50  MW/cm2) 

0.01 

0.003 

g  (I' =300  mW/cm^) 

60 

18 

g  (I'=500  mUycm^) 

100 

30 

is  an  artifact  of  earlier  AFIT  work  on  the  FRIZ  and  is  not  significant 
The  parameter,  is  assigned  a  value  several  orders  of  magnitude 

greater  than  to  reflect  the  much  lower  probability  of  a  donor 

recombination  event  versus  a  trap  recombination,  or  ionization,  event. 
As  mentioned  earlier,  some  parameters  will  be  modified  individually  in 
some  investigations. 
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Selected  results  from  the  first  set  of  validation  runs  are 


contained  in  the  Appendix,  Figures  41  to  45.  These  figures  reflect 
normalized  data  for  the  SI,  conducting,  and  DI  PRIZ  for  master  model  one 
and  parameter  set  one.  Data  is  plotted  in  a  three-dimensional  landscape 
format  versus  distance  and  time,  and  includes  density  of  ionized  donors, 
density  of  ionized  traps,  and  density  of  free  electrons.  Trends  in 
these  plots,  as  well  as  other  validation  data,  indicate  that  master 
model  one  qualitatively  duplicates  the  major  features  of  the  PRIZ,  as 
described  in  this  report  and  others.  For  example,  all  three  types  of 
PRIZ  show  identical  results  for  density  of  ionized  donors,  as  reflected 
in  Figure  41.  This  is  expected  since,  neglecting  the  relatively 
insignificant  donor  recombination  term,  the  rate  equation  that  models 
this  factor  is  unaffected  by  the  insulating  characteristics  of  the 
cathode  and  anode.  Data  for  the  density  of  ionized  traps  and  density  of 
free  electrons  clearly  shows  the  areas  of  electron  concentration  and 
enhanced  trapping  as  the  bottleneck  region  retreats  toward  the  cathode. 
The  data  also  shows  the  accelerated  movement  of  the  bottleneck  towards 
the  cathode  for  the  conducting  versus  the  SI  PRIZ,  as  can  be  seen  by 
comparisons  of  Figures  42  and  43  and  Figures  44  and  45.  This  is 
primarily  due  to  injection.  Validation  data  for  the  DI  PRIZ  showed  the 
area  of  intense  trapping  as  the  electrons  bunch  at  the  insulated  anode. 
When  the  DI  plots  were  rescaled  to  overcome  the  dominance  of  anode 
electron  concentration  and  trapping,  the  DI  PRIZ  data  for  density  of 
ionized  traps  and  density  of  free  electrons  was  indistinguishable  from 
the  corresponding  SI  PRIZ  data.  Data  for  net  charge  density  clearly 
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showed  the  formation  of  the  net  positive  space  charge  region  near  the 
cathode.  The  net  positive  space  charge  region  for  the  conducting  FRIZ 
was  obviously  smaller  in  magnitude  and  thickness  than  the  corresponding 
values  for  the  DI  and  SI  FRIZ.  As  discussed  earlier,  this  reduction  is 
due  to  injection,  and  is  especially  pronounced  at  low  intensity 
exposures.  Electric  field  data  and  total  current  data  were,  in  general, 
also  consistent  with  expected  results. 

In  addition  to  the  aspects  mentioned  above,  some  other  data  is 
presented  here  to  further  validate  master  model  one.  Figure  9,  below, 
shows  data  on  the  net  charge  density  that  accumulates  in  the  crystal 
over  time.  As  expected,  the  DI  FRIZ  accumulates  no  net  charge,  while 
the  SI  FRIZ  accumulates  a  net  positive  charge,  since  some  electrons 
escape  through  the  anode.  The  conducting  FRIZ  accumulates  a  net 
negative  charge  since,  in  this  particular  case,  more  electrons  are 
injected  through  the  cathode  than  escape  through  the  anode.  Figure  10 
shows  the  cathode  electric  field  versus  time  for  the  three  types  of 
FRIZ.  The  SI  and  DI  curves  show  that  the  cathode  field  becomes  more 
negative  as  the  positive  space  charge  accumulates  near  the  cathode.  The 
conducting  FRIZ  cathode  field  also  initially  becomes  more  negative, 
however,  the  change  is  so  small  as  to  be  imperceptible  on  the  given 
scale.  The  behavior  of  the  cathode  electric  field  for  the  conducting 
FRIZ,  at  low  intensity  exposure,  was  surprising  since  we  expected  the 
field  to  become  much  more  negative,  a  condition  that  would  encourage 
injection.  Figure  11  shows  the  photocurrent  through  the  three  types  of 
FRIZ.  As  expected,  the  conducting  FRIZ  photocurrent  is  greater  due  to 
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Figure  10.  Cathode  Electric  Field 
Master  Model  One,  Parameter  Set  Oi 


injection.  The  quick  rise  in  photocurrent  shows  that  injection  becomes 
a  factor  very  early  in  the  exposure,  as  expected  for  low  intensity 
illumination.  Figures  12  and  13  show  the  net  charge  density  versus  time 
in  bins  one  and  two,  respectively.  As  expected,  the  net  charge  density 
in  the  conducting  FRIZ  is  reduced,  versus  the  SI  and  DI  FRIZ,  due  to 
injection.  The  only  surprise  encountered  in  the  validation  was  the 
extremely  small  change  in  the  cathode  electric  field,  versus  time,  for 
the  conducting  FRIZ.  All  the  other  data,  however,  indicates  that  the 
model  is  successfully  accounting  for  the  effects  of  injection.  After 
the  validation  of  master  model  one,  computer  runs  were  made  for  master 
model  two,  parameter  set  one.  These  runs  showed  the  same  qualitative 
behavior  as  reported  for  master  model  one. 

Reproduction  of  Low  Intensity  Exposure  Data 

Next,  we  attempted  to  reproduce  the  Russian  data  reported  in 
Figures  5  and  6.  Figure  14  shows  the  results  of  our  attempts  to 
reproduce  Figure  5,  the  Russian  SI  FRIZ  data  at  low  intensity  exposure 
(I'  =  50  /iW/cm^)  ,  using  both  master  models  and  both  parameter  sets. 
Although  all  three  combinations  of  models  and  parameter  sets  produced 
virtually  the  same  results,  our  data  shows  charge  distributions  that  are 
approximately  30  times  less  than  those  reported  by  the  Russians. 
Encouragingly,  all  exposure  curves  were  qualitatively  identical  to  the 
Russians'  and  show  the  same  crossing  values  for  the  z  axis.  Also,  the 
relative  magnitudes  of  all  three  exposure  curves  are  identical  to  the 
relative  magnitudes  of  the  Russian  curves. 
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Figure  13.  Net  Charge  Density  vs.  Time,  Bin  2,  I'  =30  fiV/crn^, 
ng  =  0 .  0.1,  Master  Model  One,  Parameter  Set  One 

Figure  15  shows  the  results  of  our  attempts  to  reproduce  Figure  6, 
the  Russian  conducting  PRIZ  data  at  low  intensity  exposures.  Again,  the 
qualitative  aspects  of  the  plots  are  identical  to  the  Russians'  plots, 
and  all  three  combinations  of  master  models  and  parameter  sets  produce 
similar  data.  Crossing  values  are  the  same  as  Russian  values  as  are 
relative  magnitudes  of  the  three  exposure  curves.  Most  importantly, 
comparison  cf  Figures  14  and  15  shows  that  the  effects  of  injection  on 
the  charge  distribution  are  identical  to  effects  observed  in  the  Russian 
plots.  Specifically,  we  can  see  that  both  the  magnitude  and  thickness 
of  the  positive  space  charge  distributions  are  reduced  by  approximately 
50  percent  due  to  injection.  Despite  these  successes,  however,  our 
plots  are  once  again  off  by  a  factor  of  30.  To  improve  our  quantitative 
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Figure  14.  Net  Charge  Density  Distribution.  I'  =  50  pW/cm 
^  =  0.4  (1),  2  (2),  and  4  (3) 
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Figure  15.  Net  Charge  Density  Distribution,  I'  -  50 
r  =  0.4  (1),  2  (2),  and  4  uJ/cm^  (3) 


agreement  with  the  Russian  data,  some  of  the  more  questionable  parameter 
values  were  varied  to  determine  their  potential  effects  on  quantitative 
results.  Using  master  model  two,  parameter  set  one,  and  the  4  /iJ/cm^ 
exposure  curve  in  Figure  14  as  a  basis.  Figure  16  shows  the  effects  of 
individually  varying  values  for  absorption  coefficient  and  average 
electron  drift  time.  Figure  16  confirms  Bryksin's  assertion  in  a  1986 
publication  that,  "the  magnitude  of  the  light  absorption  coefficient  is 
decisive  as  regards  the  occurrence  of  a  particular  distribution."  (9:80) 
To  improve  our  results  by  a  factor  of  30,  however,  the  absorption 
coefficient  must  be  increased  well  beyond  the  value  shown  in  Figure  16 
(a  =  9 ,  q'  =  200  cm"^)  .  This  extreme  value  is  not  justified  for  write 
light  wavelengths  of  441  nm. 

Figure  16  also  shows  that  changing  the  value  of  average  electron 
drift  time  adversely  affects  crossing  values  without  greatly  improving 
quantitative  results.  Other  factors  investigated  were  density  of  traps, 
density  of  donors,  and  number  of  bins.  None  of  these,  however,  could 
account  for  the  discrepancies:  and  combinations  of  changes  in  some 
values,  within  reasonable  bounds,  were  also  unsuccessful.  In  spite  of 
quantitative  discrepancies,  modelling  of  all  important  qualitative 
aspects  of  the  Russian  data  was  successful  enough  to  encourage  further 
studies.  Figures  17  and  18  highlight  the  qualitative  successes  of  our 
low  exposure  intensity  modelling  efforts  by  plotting  the  Russian  data 
and  our  data  on  the  same  graphs.  Our  data  has  been  scaled  by  a  factor 
of  26  and  uses  master  model  two  with  parameter  set  one. 
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figure  16.  Effects  of  Absorption  Coefficient  and  Average  Electron  Drift 
Time  on  Net  Space  Charge  Density,  I'  -=  50  ^W/cm^,  n^  =  0 ,  ^  =4  ^J/cm^ 
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Figure  17.  Comparison  of  Russian  Low  Intensity  Exposure  Data  with  Our 
Data  Scaled  up  by  a  Factor  of  26,  I'  -  50  /iW/cm^,  ng  -=  0,  ^  -=  0.4  (1)  , 
2  (2),  and  4/iJ/cm2  (3)  (6:375) 


Figure  18.  Comparison  of  Russian  Low  Intensity  Exposure  Data  with  Our 
Data  Scaled  up  by  a  Factor  of  26,  1'  -  50  /iW/cm^,  ng  =  0.1,  ^  ■=  0.4  (1)  , 
2  (2),  and  (3)  (6:375) 


Reproduction  of  High  Intensity  Exposure  Data 

Next,  we  attempted  to  reproduce  Russian  high  intensity  exposure 
data  (I'  “  500  mW/cm^)  shown  in  Figure  7.  Landscape  plots  showed  the 
same  qualitative  features  used  to  validate  our  models  at  low  intensity 
exposure.  Data  for  net  charge  density  accumulated  in  the  crystal  versus 
time,  however,  showed  one  major  difference  from  the  low  intensity 
exposure  data.  The  net  charge  accumulated  for  a  conducting  FRIZ  was 
positive,  for  this  particular  high  intensity  exposure,  indicating  more 
electrons  were  extracted  from  the  anode  than  injected  through  the 
cathode.  This  is  consistent  with  earlier  assertions  that  injection  at 
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high  intensities  and  low  exposures  is  relatively  insignificant. 

Figure  19  shows  net  charge  accumulation  for  the  SI,  conducting,  and  DI 
FRIZ  under  high  intensity  illumination.  Figure  20  shows  the  cathode 
electric  field  versus  time  for  high  intensity  exposures.  Once  again, 
the  conducting  FRIZ  data  shows  striking  differences  from  the 
corresponding  low  intensity  exposure  data.  For  high  intensity,  the 
cathode  field  for  the  conducting  FRIZ  becomes  more  negative  at 
approximately  the  same  rate  as  that  for  the  SI  and  DI  FRIZ.  This 
indicates  that  injection  has  not  yet  become  significant  enough  to  modify 
evolution  of  the  space  charge,  as  happens  almost  immediately  for  low 
intensity  exposures.  Figure  21  shows  that  the  photocurrent  through  the 
FRIZ  at  high  intensity  exposure  is  identical  for  the  SI  and  conducting 
FRIZ.  Likewise,  Figures  22  and  23  show  that  the  net  charge  density 
versus  time  in  bins  one  and  two  are  almost  identical.  All  of  these 
features  indicate  that  injection  is  relatively  insignificant  at  high 
intensity  exposure. 

Figure  24  shows  our  attempts  to  reproduce  Russian  data  reflected  in 
Figure  7.  Once  again,  the  qualitative  features  of  the  plot  are  similar 
to  the  Russian  data.  In  this  case,  however,  our  magnitudes  are  only 
approximately  seven  times  less  than  the  Russians'.  Another  significant 
difference  is  the  fact  that  parameter  set  two  produces  significantly 
different  data  than  parameter  set  one.  The  z  axis  crossings  for 
parameter  set  two  are  closer  to  the  Russian  values  than  crossings  for 
parameter  set  one,  although  neither  set  produces  crossings  as  close  to 
the  Russian  values  as  we  experienced  in  reproducing  low  intensity 
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Figu:  3  19.  Net  Charge  Accumulated  vs.  Time,  I'  =  500  mW/cm 
0.1,  “'.aster  Model  One,  Parameter  Set  One 
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Figure  20.  Cathode  Electric  Field  vs.  Time,  I'  -  500  mW/cm' 
0.1,  Master  Model  One,  Parameter  Set  One 
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Figure  21.  Photocurrent  vs.  Time,  I'  =  500  mW/cm^,  ng  -  0 ,  0.1, 
Master  Model  One,  Parameter  Set  One 
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Figure  22.  Net  Charge  Density  vs.  Time,  Bin  1,  I'  =  500  mW/cm^, 
Hg  =  0 ,  0.1,  Master  Model  One,  Parameter  Set  One 


Figure  23.  Net  Charge  Density  vs.  Time,  Bin  2,  I'  -=  500  mW/cm^, 
ns  =  0 ,  0.1,  Master  Model  One,  Parameter  Set  One 


exposure  data.  The  shift  of  crossings  to  higher  z  values  for  parameter 
set  two  is  primarily  due  to  the  increased  electron  mobility 

=  0.03  cm^/V*sec  for  set  two,  ^  -0,01  cm^/V»sec  for  set  one),  which 
allows  the  photogenerated  electrons  to  drift  more  rapidly,  enhancing  net 
positive  charge  development.  Figure  25  shows  plots  of  the  Russians'  and 
our  high  intensity  exposure  data  on  the  same  graph.  Our  data  has  been 
scaled  upward  by  a  factor  of  seven  and  reflects  data  for  a  conducting 
PRIZ  using  master  model  two  and  parameter  set  one.  Again,  the  reasons 
for  the  discrepancies  in  quantitative  results  are  unclear.  The  fact 
that  the  discrepancy  is  seven  for  high  intensity  plots  and  30  for  low 
intensity  plots  is  inexplicable.  As  for  low  intensity  exposures. 


65 


I 


-20  - - - — - - - 1 - 1 - - 

0.00  0.05  O.IO  0.15  0.20  0.25  0.30 


Figure  24.  Net  Charge  Density  Distribution,  1'  -  500  mW/cm^,  n^  -  0.1, 
r  =  4  (1),  8  (2),  12  (3),  24  pJ/cm^  (4) 
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Figure  25.  Comparison  of  Russian  High  Intensity  Exposure  Data  with  Our 
Data  Scaled  up  by  a  Factor  of  7,  I'  ■=  500  mW/cm^,  ng  =  0.1,  ^  =  4  (1), 

8  (2),  12  (3),  and  24  nJ /cm^  (4)  (6;375) 


attempts  were  made  to  modify  parameters,  within  reasonable  bounds,  to 
account  for  the  factor  of  seven  discrepancy.  Plots  with  individually 
varied  parameters  such  as  number  of  bins,  absorption  coefficient, 
average  electron  drift  time,  density  of  donors,  and  density  of  traps, 
however,  could  not  account  for  the  discrepancies. 

Next,  we  made  attempts  to  demonstrate  the  1988  paper's  assertions 
that  data  in  Figure  7  differs  by  less  than  one  percent  for  an  SI  PRIZ 
(6:375).  Since  we  could  not  quantitatively  reproduce  their  data  in 
Figure  7,  we  had  to  rely  on  comparisons  for  an  SI  and  a  conducting  PRIZ 
using  data  we  had  generated  in  Figure  24.  Figure  26  demonstrates  that 
the  assertion  is  accurate,  using  our  data  from  master  model  two  and 
parameter  set  one.  Though  data  for  a  DI  PRIZ  showed  significant 
deviations,  the  SI  and  conducting  PRIZ  data  for  all  exposures  shown  is 
virtually  identical.  Figure  26  reinforces  previous  assertions 
concerning  the  relative  insignificance  of  injection  at  high  intensities 
and  low  exposures. 

Finally,  we  analyzed  the  1988  paper's  assertion  that,  "As  the 
intensity  of  the  writing  light  increases,  the  range  of  exposures  at 
which  injection  has  only  insignificant  effect  becomes  broader."  (6:375) 
Figures  27  to  29  show  plots  of  net  charge  density  versus  z  for  exposures 
of  4  and  40  ;^J/cm^  and  intensities  of  30  /iW/cm^,  3  mW/cm^,  and 
300  mW/cm^,  respectively.  These  plots  show  that  the  assertion  given 
above  is  also  accurate,  using  data  from  master  model  two  and  parameter 
set  one. 
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Figure  26.  Comparison  of  High  Intensity  Exposure  Data  for  an  SI  and  a 
Conducting  PRIZ,  I'  -  500  mW/cm^,  ng  -  0,  0.1,  ^  -  4  (1)  ,  8  (2)  ,  12  (3), 
and  24  ^ij/cm^  (4) 
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Figure  27.  Comparison  of  Low  Intensity  Exposure  Data  for  an  SI  and  a 
Conducting  PRIZ,  I'  =  30  ^W/cm^,  ng-O,  0.1,  ^  >=4  and  40  /iJ/cm^ 
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Figure  28.  Comparison  of  Medium  Intensity  Exposure  Data  for  an  SI  and 
a  Conducting  FRIZ,  I'  =  3  mW/cm^,  ng^O,  0.1,  ^  -  4  and  40  pJ/cm^ 
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Figure  29.  Comparison  cf  High  Intensity  Exposure  Data  for  an  SI  and  a 
Conducting  FRIZ,  I'  =  300  mW/c'^^,  ng  -  0,  0.1,  ^  -  4  and  40 
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Investigations  of  Reciprocity 


Before  attempting  to  reproduce  reciprocity  data,  we  analyzed  the 
development  of  the  net  charge  distribution,  of  the  SI  and  conducting 
PRIZ,  after  turn-off  of  the  write  light.  Data  for  the  DI  PRIZ  was 
virtually  identical  to  that  for  the  SI  PRIZ  and  will  not  be  discussed 
separately.  When  the  write  light  is  turned  off,  the  evolution  of  the 
net  charge  distribution  depends  on  the  injecting  characteristics  of  the 
cathode,  material  properties  of  the  PRIZ  crystal,  such  as  electron 
mobility,  and  lag  time  effects.  Lag  time  effects  are,  as  previously 
discussed,  more  prevalent  at  high  intensities  and  low  exposure  times, 
where  the  charge  distribution  has  not  had  time  to  fully  develop. 

Analysis  was  accomplished  by  continuing  the  calculations  of  the 
differential  equations  after  setting  the  generation  coefficient  to  zero. 
Results  for  all  low  intensity  runs  (I'  =  30  ^W/cm^)  showed  that  the  net 
space  charge  distributions  did  not  change  perceptibly  after  the  PRIZ  was 
allowed  to  decay,  or  relax,  for  up  to  10  time  constants  (10  •  tg'  =  1 
msec)  past  write  light  turn-off.  This  held  true  for  the  SI,  conducting, 
and  DI  PRIZ  and  reflects  the  inherent  optical  storage  capabilities  of 
the  PRIZ.  Basically,  for  low  intensity  illumination,  the  characteristic 
exposure  times  allow  the  charge  separation  to  continually  develop 
throughout  the  exposure  process.  Therefore,  when  the  write  light  is 
turned  off,  little  additional  charge  separation  is  left  to  occur.  For 
high  intensity  illumination,  however,  the  relaxation  process  was  more 
evident,  as  illustrated  in  Figures  30  and  31  for  the  SI  and  conducting 
PRIZ,  respectively.  The  net  space  charge  distribution  in  the  SI  PRIZ 
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widens  at  a  progressively  slower  rate  after  turn-off  of  the  write  light. 
Most  of  the  additional  net  positive  space  charge  distribution,  however, 
occurs  during  the  first  time  constant  past  write  light  turn-off.  Our 
analysis  only  considered  relaxation  times  up  to  ten  time  constants. 
Eventually,  of  course,  the  positive  space  charge  distribution  will  decay 
due  to  thermal  ionization,  which  is  represented  in  master  model  one,  and 
other  equilibrium  processes.  Relaxation  in  the  conducting  FRIZ  is 
significantly  different  as  evident  from  Figure  31.  After  turn-off  of 
the  write  light  the  net  positive  space  charge  region  initially  expands 
hut  quickly  contracts  again.  The  maximum  net  space  charge  region 
appears  to  occur  after  one  time  constant. 

Referring  back  to  Figure  8,  the  relaxation  characteristics 
mentioned  above  could  explain  why  the  curves  for  high  and  low  intensity 
exposures  do  not  exhibit  reciprocity,  using  our  criterion  for  total  area 
under  the  net  positive  space  charge  curve.  If  response  measurements  are 
made  one  time  constant  following  the  completion  of  a  given  exposure, 
then  the  low  intensity  curve  will  be  unaffected  while  the  high  intensity 
curve  will  expand  out  to  encompass  more  total  area.  In  doing  so,  the 
resulting  curves  could  meet  our  reciprocity  criterion.  Unfortunately, 
however,  reproduction  of  Figure  8  using  our  data  for  high  and  low 
intensif.'  expcc’uies  did  not  even  qualitatively  resemble  the  Russian 
data,  as  snown  'n  Figure  32.  Our  data  indicates  that  the  curves  exhibit 
reciprocity  without  accounting  for  relaxation.  iieciprocity  fails,  in 
our  data,  if  a  ore  time  constant  delay  is  used  prior  to  measuring 
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Figure  30.  Relaxation  of  Net  Positive  Space  Charge  Following  Turn-off 
of  Write  Light,  1'  =  500  mW/cm^,  ng  =  0,  ^  *=24  /ij/cm^ 
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Figure  31.  Relaxation  of  Net  Positive  Space  Charge  Following  Turn-off 
of  Write  Light,  I'  =  500  mW/cm^,  nj  =  0.1,  ^  -  24  /jJ/cm^ 
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Figure  32.  Relaxation  Distribution  of  Net  Space  Charge  for  High 
(500  mW/cm^)  and  Low  (50  ^W/cm^)  Intensity  Exposures,  ng  =  0.1, 

^  =4 


response  Fii'r>-hcr  analysis  indicates  that  the  apparent  reciprocity 
shown  by  our  data  in  Figure  32  is  not  indicative  of  an  overall  trend. 

To  further  investigate  reciprocity,  we  analyzed  net  space  charge 
density  data  for  an  SI  and  a  conducting  FRIZ  at  write  light  intensities 
from  30  ^W/cm^  to  300  mW/cm^  and  exposures  from  4  to  40  ;iJ/cm^.  The 
intensities  correspond  to  those  in  the  1988  paper  for  which  reciprocity 
was  experimentally  analyzed,  and  the  exposures  correspond  to  values 


within  and  without  the  experimentally  established  region  of  reciprocity. 
Net  charge  density  data  was  recorded  for  both  types  of  PRIZ,  and  for 
several  intensities  and  exposures,  at  times  corresponding  to  zero  delay 
and  one  time  constant  delay.  The  results  for  the  SI  PRIZ,  using  no 
delay  in  recording  the  net  space  charge  density,  are  shown  in  Figure  33. 
This  figure  indicates  that  the  SI  PRIZ  with  no  delay  is  not  reciprocal 
at  the  three  lower  exposure  values  (4,  12,  and  20  but  becomes 

reciprocal  at  the  higher  exposure  value  of  40  fiJ/cm^.  Figure  34, 
however,  indicates  that  the  SI  PRIZ  exhibits  reciprocity  at  all  exposure 
and  intensity  values  when  a  one  time  constant  delay  is  used.  Figure  35 
shows  that  the  conducting  PRIZ,  with  no  time  delay,  thoroughly  fails  our 
reciprocity  criterion.  Figure  36  indicates  that  the  conducting  PRIZ 
with  a  one  time  constant  delay  fails  even  more  thoroughly. 

Reciprocity  appears  to  defy  description  in  one  dimension.  One 
possible  explanation  for  our  failures  in  this  respect  is  contained  in  a 
1984  publication  by  Bryksin.  Bryksin  states,  "We  emphasize  that  the 
negative  space  charge  region  has  a  greater  influence  than  the  positive 
space  charge  region  on  the  integrated  characteristic  for  the  transverse 
field  because  it  is  located  farther  away  from  the  negative  metal 
electrode."  (12:881)  On  the  other  hand,  Petrov  states  in  a  1989  paper 
that  "the  total  positive  charge  is  usually  considerably  larger  than  the 
negative  one,  and  thus  we  can  justifiably  neglect  the  readout  light 
modulation  via  the  negative  charge  field."  (32:337)  Although  these  two 
statements  appear  contradictory,  the  specific  contexts  in  which  they  are 
applied  remove  the  ambiguity.  Bryksin's  statement  applies  to  space 
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Figure  36.  Distribution  of  Net  Space  Charge  in  a  Conducting  FRIZ  with 
a  One  Time  Constant  Delay  in  Recording  Response,  nc  -  0.1,  ^  —  4  (a) 
12  (b) ,  20  (c),  and  40  nJ/cm^  (d) 


charge  distributions  that  have  significant  regions  of  both  positive  and 
negative  space  charge.  Figure  18  on  page  61  shows  such  a  situation  for 
low  intensity  exposure.  Petrov  is  referring  to  distributions  in  which 
the  positive  space  charge  region  is  dominant.  Figure  25  on  page  66 
shows  this  situation  for  high  intensity  exposure.  The  specific  nature 
of  the  space  charge  distribution  determines  the  approximations  we  can 
use  to  predict  the  response.  Our  simple  criterion  for  predicting 
response,  based  on  the  total  area  of  the  positive  space  charge  region, 
is  inadequate  for  distributions  with  a  significant  negative  space  charge 
region . 

Investigations  of  Self  Erasure 

The  first  figures  in  both  the  1987  and  1988  papers  indicate  a 
phenomenon  called  self  erasure.  Self  erasure  occurs  when  the  PRIZ 
response,  due  to  a  stationary  image,  develops  peaks  and  then  starts  to 
diminish.  In  some  cases,  the  image  is  completely  erased  despite 
continuous  imaging.  Figure  3  depicts  the  phenomenon  of  self  erasure  in 
a  conducting  PRIZ.  Figure  4  shows  total  exposure  values  (^max) • 
various  intensities,  at  which  the  response  of  the  PRIZ  peaks.  Exposures 
beyond  these  ?n,ax  points  demonstrate  the  onset  of  self  erasure.  The 
exposure  times  to  produce  self  erasure  are  much  greater  than  any  we  have 
analyzed  so  far.  For  example,  at  50  ^W/cm^,  the  greatest  exposure  time 
we  have  analyzed  up  until  now  is  80  msec  ~  4  /jJ/cm^)  .  To  experience 

self  erasure  at  this  intensity,  we  must  use  exposure  times  greater  than 
one  second.  For  500  mW/cm^,  we  have  analyzed  exposure  times  of  48  ^sec 
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=  Ik  pj/cm^).  To  experience  self  erasure  at  this  intensity,  we  must 
use  exposure  times  greater  than  one  millisecond. 

To  analyze  self  erasure  in  one  dimension,  we  used  intensities  of 
30  ^W/cm^  and  300  mW/cm^  with  exposure  times  of  2.7  sec  =  80  /iJ/cm^) 
and  16,2  msec  =4.9  mJ/cm^) ,  respectively.  These  exposure  values 

should  reach  well  into  the  self  erasure  region  according  to  data  in 
Figure  4.  Once  again,  to  model  the  two-dimensional  transverse  electric 
field  effects  using  our  one-dimensional  models,  we  intuitively  analyzed 
the  fields  produced  when  our  one-dimensional  net  space  charge 
distribution  is  placed  adjacent  to  a  net  zero  space  charge  region.  As 
we  found  in  our  discussions  on  reciprocity,  the  overall  response 
produced  by  this  configuration  is  very  complicated.  However,  if  we  can 
show  that  our  charge  distribution  approaches  zero  as  the  exposure 
increases,  then  the  transverse  fields  between  the  two  adjacent  regions 
would  also  approach  zero.  The  absence  of  transverse  fields  would  result 
in  zero  response,  indicating  self  erasure.  Certainly,  a  uniformly  zero 
charge  distribution  is  not  the  only  way  to  achieve  zero  response.  A 
favorable  combination  of  positive  and  negative  space  charge  regions 
throughout,  the  length  of  the  crystal  could  conceivably  produce 
transverse  fields  that  result  in  zero  net  response.  This  situation, 
however,  is  much  more  difficult  to  analyze  with  one-dimensional  models. 

Intuitively,  we  anticipated  that,  due  to  the  bottleneck  effect,  the 
net  positive  space  charge  region  would  retreat  toward  the  cathode  as  the 
exposure  increased.  Whether  the  bins  nearest  the  cathode  would  actually 
diminish  to  zero  net  charge  density,  however,  was  unclear.  Figure  37 
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indicates  that  all  bins  except  the  first  one  do  indeed  go  to  zero  as  the 
exposure  increases.  Unfortunately,  the  net  charge  density  in  bin  two 
passes  through  the  zero  value  and  becomes  increasingly  negative.  The 
response  produced  by  such  a  distribution  is  difficult  to  predict. 

Figure  37  was  produced  using  master  model  one,  parameter  set  one,  and  an 
intensity  of  300  mW/cm^.  Since  the  1988  paper  never  specified  material 
parameters  for  their  experiment,  however,  we  varied  several  of  the 
values  in  parameter  set  one  in  attempts  to  drive  all  bins  to  zero  net 
charge  density.  Figure  38  shows  the  results  of  increasing  the  cathode 
injection  coefficient,  ng ,  from  0.1  to  10.  As  can  be  seen,  the  charge 
in  bin  two  is  significantly  reduced,  although  it  still  eventually 
becomes  negative.  Decreasing  the  average  electron  drift  time,  r^,  to 
0.009  versus  0.9  produced  data  identical  to  that  for  ng  =  10.  Although 
we  were  able  to  affect  the  charge  density  in  bin  two  by  adjusting 
parameter  values,  we  were  never  able  to  arrest  the  linear  rise  of  the 
charge  density  in  bin  one.  Interestingly,  individual  plots  of  density 
of  ionized  trap  data  from  master  model  one  indicated  that  filling  of 
traps  occurred  throughout  the  length  of  the  crystal  during  the  course  of 
the  exposure.  When  we  switched  to  master  model  two,  which  does  not 
model  the  depletion  of  traps,  however,  we  saw  little  change  in  the  net 
charge  densities  in  bins  one  and  two. 

Figure  39  shows  data  for  master  model  one,  parameter  set  one,  and 
an  intensity  of  30  piW/cm^ .  Here,  the  net  charge  density  in  bin  one  did 
not  continually  increase  at  a  linear  rate,  although  it  also  never 
decreased.  The  charge  in  bin  two  never  goes  to  zero  and  appears  to 
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Net  Charge  Density  vs.  Time  in  Bins  1  (1),  2  (2),  3  (3), 
'  =  300  mW/cm^,  ns“0.1,  2”  -  0  to  4.9  mJ/cm^ 


Net  Charge  Density  vs.  Time  in  Bins  1  (1),  2  (2),  3  (3), 
'  =  300  mW/cm^,  ng  -  10,  ^  -  0  to  4.9  mJ/cm^ 
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level  out  at  a  moderately  high  value  of  net  positive  charge  density. 

The  charge  densities  in  higher  bins  approach  and  sometimes  cross  zero, 
then  remain  at  relatively  small  positive  or  negative  values.  Figure  40 
shows  that  by  switching  to  master  model  two,  we  were  able  to  force  the 
value  of  net  charge  density  in  all  bins,  except  bin  one,  towards  zero. 
More  significantly,  we  were  also  able  to  cause  the  net  charge  density  in 
bin  one  to  level  off,  although  it  still  never  decreased.  The 
differences  produced  by  the  two  master  models  indicated  that 
recombination  of  donors,  filling  of  traps,  or  detrapping  was  having  a 
major  affect  on  the  results  from  master  model  one.  None  of  these 
effects  was  observed,  however,  from  individual  plots  of  density  of 
ionized  donors  or  density  of  ionized  traps.  More  perplexing  was  the 
fact  that  the  data  showed  no  sensitivity  to  changes  in  nj  from  0.01  to 
10.  Plots  of  the  cathode  electric  field  versus  time  revealed  that,  once 
again,  the  cathode  field  was  staying  at  or  near  its  initial  value,  a 
condition  that  does  not  allow  significant  injection.  This  recurring 
feature  of  our  analytical  data,  for  lew  intensity  exposures,  prompted  a 
more  detailed  analysis  of  the  specific  conditions  for  which  the  feature 
was  observed.  This  analysis  showed  that  when  the  ratio  of  the  injection 
coefficient,  ng ,  to  a  quantity  that  characterizes  the  rate  of  charge 
separation,  gr^,  is  greater  than  one,  the  cathode  field  exhibited  little 
variation  from  its  initial  value.  Bryksin  previously  characterized  the 
behavior  of  his  analytical  models  with  respect  to  the  parameter  ng/gr^j 
in  several  papers.  In  his  earlier  models,  he  established  that  values 
for  this  parameter  must  be  less  than  one  for  accurate  results  (16:1350). 


83 


Figure  39. 
and  A  (4)  , 
One 


^  0  5000  10000  15000  20000  25000  30000 

t 
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Figure  40.  Net  Charge  Density  vs.  Time  in  Bins  1  (1),  2  (2),  3  (3), 
and  4  (4),  I'  =  30  /jW/cm^,  ng  =  0.1,  ^  -  0  to  80  /jJ/cm^,  Master  Model 


Though  subsequent  models  removed  the  assumption  that  mandated  this 
restriction,  the  behavior  of  the  FRIZ  is  still  tied  to  the  value  of  this 
ratio  relative  to  one.  For  both  parameter  sets,  this  means  that  a 
generation  coefficient  (g)  value  of  approximately  0.1  corresponds  to  the 
borderline  region  where  ng/gr^j  =  1.  This  borderline  region  correlates 
well  with  our  results,  and  should  account  for  the  observed  differences 
in  FRIZ  dynamics  between  low  intensity  exposures  (ng/gr^j  =  10)  and  high 
intensity  exposures  (ng/gr^j  ~  0.001).  Therefore,  the  parameter 
ng/gTd  =  1  defines  a  critical  transition  point  in  FRIZ  dynamics ,  and 
must  be  accounted  for  in  any  interpretation  of  analytical  results  from 
our  models.  Accepting  this,  as  well  as  the  basic  validity  of  our  model, 
we  can  easily  see  why  the  cathode  electric  field  never  varies 
significantly  from  its  initial  value  for  low  intensity  exposures. 

Though  photoionization,  electron  drift,  and  electron  trapping  enjoy  an 
extremely  short  injection-free  period  in  which  to  build  a  cathode 
electric  field,  the  injection  current  quickly  becomes  dominant  and  then 
relaxes.  The  reasons  the  cathode  field  subsequently  becomes  less 
negative,  indicating  negative  injection,  are  unclear.  An  adjustment  of 
the  boundary  conditions  for  bin  one  or  a  modification  of  the  equation 
modelling  the  injection  current  might  eliminate  this  negative  injection. 

In  spite  of  the  difficulties  mentioned  above,  we  were  still  able  to 
show  that  the  net  charge  densities  in  all  bins,  except  bin  one,  approach 
zero  as  the  exposure  continues.  An  intuitive  analysis  of  the  charge  in 
bin  one  indicates  that  the  only  ways  for  this  last  remnant  of  net 
positive  space  charge  to  diminish  is  through  donor  recombination  or 
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enhanced  trapping.  Both  of  these  events  are  dependent  on  the  density  of 

free  electrons  in  bin  one.  Analysis  of  the  data  for  the  normalized 

density  of  electrons  in  bin  one,  however,  shows  a  minuscule  value  of 

0.1,  that  is  achieved  almost  immediately,  and  does  not  appear  to  vary  at 

all  over  time.  In  comparison,  the  normalized  density  of  electrons  in 

bin  2  climbs  to  8000  over  the  course  of  the  exposure.  Although  we 

appear  unable  to  affect  this  last  remnant  of  charge  in  bin  one,  the 

following  discussion  indicates  that  we  might  be  able  to  simply  ignore 

it,  and  thereby  claim  success  in  one-dimensionally  demonstrating  self 

erasure.  In  a  1984  paper  by  Bryksin,  he  states,  "The  contribution  from 

the  positive  space  charge  region  to  the  transverse  field  for  small  z  is 

therefore  shielded  more  strongly,  owing  to  mirror  image  forces  on  the 

negative  electrode."  (12:881)  Petrov  elaborates  on  this  shielding 

effect  near  the  cathode  as  follows: 

Contraction  of  the  charge  layer  towards  the  electrode 
accompanied  by  slowing  down  of  the  charge  density  growth 
causes  a  decrease  of  the  readout  light  modulation  amplitude  by 
virtue  of  the  fact  that  the  transverse  components  of  the 
charge  field  near  the  electrode  turn  out  to  be  weakened  to  a 
great  extent  by  the  field  of  the  "mirror"  charge  produced  by 
the  electrode.  (32:347) 

If  we  accept  these  statements,  then  the  charge  density  in  bin  one,  and 
perhaps  in  bin  two,  should  not  modulate  the  read  light  to  any  great 
degree.  Therefore,  the  response  should  eventually  diminish  with 
increasing  exposure,  indicating  self  erasure.  Similar  studies,  using  a 
progressively  greater  number  of  bins,  could  be  used  to  confirm  that  the 
net  positive  space  charge  always  retreats  into  the  progressively  smaller 
first  bin.  In  the  limit  of  infinitesimal  bin  size,  a  residual  net 


86 


positive  space  charge  exclusively  in  bin  one  would  most  definitely 
indicate  self  erasure  in  one  dimension.  Our  analysis  of  self  erasure 
neglects  the  catastrophic  effects  of  dielectric  breakdown  as  reported  by 
Anderson  (1:1).  For  high  intensity  self  erasure  data,  we  are  well  past 
the  exposure  value  ~  1.0  mJ/cm^)  where  this  damage  occurs  (23:21). 

We  did  not  analyze  the  SI  PRIZ  at  extremely  high  exposure  values. 
Such  an  exercise  could  help  isolate  the  effects  of  injection. 
Intuitively,  we  would  expect  the  net  positive  space  charge  region  in  the 
SI  PRIZ  to  also  retreat  to  the  first  few  bins,  albeit  at  a  slower  rate. 
If  these  bins  can  indeed  be  disregarded,  then  the  SI  PRIZ,  with  no 
injection,  would  also  exhibit  self  erasure. 
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VI .  Conclusions  and  Recommendations 


Our  attempts  to  analytically  reproduce  qualitative  aspects  of 
Russian  analytical  data  in  the  1988  paper  were  highly  successful.  The 
quantitative  differences  in  our  analytical  data  and  Russian  analytical 
data,  however,  are  inexplicable.  Our  models  appear  to  accurately 
reproduce  the  effects  of  injection  in  the  conducting  FRIZ,  however,  some 
data  was  troubling.  For  example,  the  cathode  electric  fields  at  low 
intensity  exposures  consistently  went  less  negative,  indicating  negative 
injection.  This  occurred  in  spite  of  the  fact  that  the  corresponding 
net  space  charge  distribution  indicated  conditions  favorable  for  the 
formation  of  an  increased  cathode  electric  field  which  should,  in  turn, 
result  in  injection.  Other  data  indicated,  however,  that  injection  was 
indeed  happening  in  spite  of  this  apparent  discrepancy.  The  extremely 
low  values  for  free  electron  density  in  bin  one,  at  extremely  high 
exposure  values,  were  also  surprising.  Overall,  however,  our  models 
appear  to  require  little  modification.  If  there  are  errors,  they 
probably  reside  in  our  modelling  of  bin  one  or  its  boundary  conditions. 

Reproduction  of  Russian  experimental  data  was  much  less  successful. 
Since  experimental  data  on  reciprocity  and  self  erasure  emerged  as 
essentially  multi-dimensional  effects,  we  were  not  surprised  that  our 
one-dimensional  models  were  somewhat  inadequate.  The  Russians'  failure 
to  report  specific  experimental  data  and  analytical  parameters  also 
hampered  our  efforts.  Great  benefit  could  be  derived  by  repeating  the 
Russian  experiments,  reported  in  the  1987  and  1988  papers,  using  FRIZ 
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crystals  with  known  material  parameters.  Analysis  of  the  results,  using 
the  models  specified  in  this  paper,  might  then  prove  more  conclusive. 
Future  researchers  could  also  refine  our  models  by  including  diffusion 
current,  in  addition  to  conduction  and  displacement  current,  in  the  term 
for  total  current.  Refinements  in  injection  modelling  could  resolve 
apparent  inconsistencies  in  the  behavior  of  the  injection  current  for 
low  intensity  exposures.  This  might  also  result  in  more  believable 
demonstrations  of  reciprocity  and  self  erasure.  A  more  sophisticated 
criterion  for  predicting  transverse  electric  fields,  and  subsequent 
response,  would  also  help  in  analyzing  reciprocity  and  self  erasure. 

Our  crude  criterion  is  rendered  virtually  useless,  in  many  situations, 
because  of  its  total  disregard  of  the  negative  space  charge  region  and 
the  effects  of  the  image  charge. 

The  most  obvious  lesson  learned  from  this  report  is  that  one¬ 
dimensional  modelling  is  quickly  approaching  the  point  of  diminishing 
returns.  The  very  nature  of  the  transverse  electrooptic  effect  mandates 
the  development  of  more  sophisticated  two-dimensional  models.  Accepting 
this,  the  computing  power  required  by  future  AFIT  FRIZ  resear :'^iers  can 
be  expected  to  rise  dramatically.  Software  upgrades  will  also  be 
required  since  even  our  one-dimensional  models  often  taxed  the  abilities 
of  existing  software. 
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Appendix:  Landscape  Plots 
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Figure  45.  Density  of  Free  Electrons  vs.  Time  and  Distance, 
Conducting  PRIZ,  I'  -  50  /iW/cm^,  ng  —  O.!,  ^  -  0to4  /iJ/cm^ 
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